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We derive a new version of the non-Abelian Stokes theorem for the Wilson loop in 
the SU{N) case by making use of the coherent state representation on the coset space 
SU{N)/U{1)^-'^ = Fn-1, the flag space. We consider the SU{N) Yang-Mills theory in the 
maximal Abelian gauge in which SU{N) is broken down to (7(1)^~^ . First, we show that the 
. Abelian dominance in the string tension follows from this theorem and the Abelian-projected 

effective gauge theory that was derived by one of the authors. Next (but independently), 
' combining the non-Abelian Stokes theorem with a novel reformulation of the Yang-Mills the- 

^ , ory recently proposed by one of the authors, we proceed to derive the area law of the Wilson 

■ loop in four-dimensional SU{N) Yang-Mills theory in the maximal Abelian gauge. Owing to 

' dimensional reduction, the planar Wilson loop at least for the fundamental representation 

I in four-dimensional SU{N) Yang-Mills theory can be estimated by the diagonal (Abelian) 

Wilson loop defined in the two-dimensional CP'^~^ model. This derivation shows that the 
fundamental quarks are confined by a single species of magnetic monopole. The origin of 
' the area law is related to the geometric phase of the Wilczek-Zee holonomy for U{N ~ 1). 

, The calculations are performed using the instanton calculus (in the dilute instanton-gas 

approximation) and using the large A'^ expansion (in the leading order). 



Qh! ^1. Introduction 

^ ■ 

Gell-Mann and Zweig& predicted in the mid-1960s tliat all hadrons (i.e., baryons 
and mesons) are composed of the fundamental constituents having fractional charges, 
^ ' "l"!^' with e being the elementary charge. Now, the fundamental constituent 

. is called the quark and the proposed theory is called the quark model. The predic- 

tions of this model are consistent with the results of experiments performed over the 
past thirty years. The strong interaction among quarks and anti-quarks is mediated 
by the gluon, which is described by the quantized Yang-Mills gauge field theory. 0) 
The present fundamental theory describing the quark and the gluon is provided by 
quantum chromo dynamics (QCD), which is a non-Abelian gauge theory with the 
gauge group G = SU{3) corresponding to three colors. However, neither an iso- 
lated quark nor an isolated anti-quark has ever been observed experimentally. In 
the present understanding, they are believed to be confined in hadrons. This is the 
hypothesis of quark confinement. Quark confinement could be explained theoreti- 
cally within the framework of QCD, although no one has achieved a rigorous proof 
of quark confinement. This is one of the most important problems to be solved in 
theoretical physics. 

QCD has a remarkable property, called asymptotic freedom, which was discov- 
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ered by Gross, Wilczek and Politzer and independently by 't Hooft.0^ Asymptotic 
freedom does not appear in the most successful quantized field theory, quantum 
electrodynamics (QED). As is well known, QED is the Abelian gauge theory for the 
electron and the photon in which the electromagnetic interaction is described by 
the quantized Maxwell gauge field theory with the gauge group G = U{1). Asymp- 
totic freedom is a consequence of gluon self-interactions. Therefore, this is a very 
characteristic feature of non-Abelian gauge theory. 

The purpose of this article is to demonstrate quark confinement within QCD 
based on the Wilson criterion for quark confinement, i.e., the area law of the 
Wilson loop. The Wilson loop is a gauge invariant quantity and hence the Wilson 
criterion is also gauge invariant. The formulation of lattice gauge theory proposed 
by Wilson is manifestly gauge invariant and does not need the gauge fixing. It is 
easy to show that the strong coupling expansion in lattice gauge theory leads to the 
area law of the Wilson loop. However, this result has not yet been continued to the 
weak coupling region, where the string tension is expected to obey the scaling law 
suggested from the result of the renormalization group based on loop calculations. 
The first indication of the area law of the Wilson loop for arbitrary coupling constant 
was foundin a study based on numerical simulations within the lattice gauge theory 
by Creutzu) for G = SU{2) and SU{3). Although the numerical evidence of quark 
confinement was indeed a great progress toward a complete understanding of quark 
confinement, the analytical proof is still lacking. 

This work was initiated to justify the dual superconductor picture of the QCD 
vacuum proposed in the mid-1970sB^ within the framework of continuum quantum 
field theory. For dual superconductivity to occur, magnetic monopoles must be 
condensed, just as ordinary superconductivity requires condensation of Cooper pairs. 
In fact, the importance and the validity of taking into account magnetic monopoles 
in quark confinement has been demonstrated, at least for simplified four-dimensional 
and lower-dimensional models, especially, by Polyakov,B and recently for the four- 
dimensional |Yang- Mills theory and QCD with extended supersymmetries by Seiberg 
and Witten.B^ In this scenario, quark confinement is realized due to condensation of 
magnetic monopoles. Recent developments in numerical simulations on the lattice!-* 
confirm the existence of dual superconductivity in QCD, at least, under a specific 
gauge fixing called the Abelian gauge. HHP 

This article gives a detailed exposition of the results on quark confinement that 
were announced in a previous articlelllP for G = SU{3) together with new results 
for G = SU{N) with arbitrary A^. They are extensions of the analyses of the 
Ya ng- Mills theory in the maximal Abelian (MA) gauge given in a series of articles, 



where the case of SU{2) was explicitly worked out. In this 
process, we have found that the extension from SU{2) to SU{3) is non-trivial, but 
the extension from SU{3) to SU{N), > 3, is rather straightforward. New features 
come out when we begin to analyze the SU (N) case with A^ > 3. It seems that 
they have been overlooked to this time in the conventional approach based on the 
maximal Abelian gauge. 

The MA gauge is a partial gauge fixing from the original non-Abelian gauge 
group G to its subgroup in which the gauge degrees of freedom of the coset 
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G/H are fixed. Even after the MA gauge, tliere is a residual gauge group H whicli 
is taken to be the maximal torus subgroup H = U{\)^~^. After the MA gauge, the 
magnetic monopole is expected to appear, since the Homotopy group it2{G/H) is 
non-trivial, i.e., 

7r2(5[/(A^)/C/(l)^-i) = ^i(C/(l)^-i) ^ 2^-1^ (1^1) 

This implies that the breaking of gauge group G ^ H hy partial gauge fixing leads 
to {N — 1) species of magnetic monopoles. However, we do not necessarily need 
to consider the maximal breaking SU{N) — > C/(l)^~^, although the maximal torus 
group is desirable as a gauge group of the low-energy effective Abelian gauge theory. 
i3> Actually, even if we restrict H to a continuous subgroupQ of G, there are other 
possibilities for choosing H, e.g., we can choose a subgroup H such that 

G D H D H ■.= U{1)^~\ (1.2) 

The possible number of cases for choosing H increases as N increases. We have 
found tiP that the group H may depend on the representation to which the quark 
belongs when > 3 and that it suffices to take H = U{N — 1) for the fundamental 
quark to be confined in the sense of the area law of the Wilson loop under the 
partial gauge fixing. Here H is equal to the maximal stability group specified by the 
highest-weight state of the representation of the quark in the Wilson loop. This is 
a new feature which does not show up in the SU{2) case. Nevertheless, this does 
not mean that the choice of the maximal torus does not lead to quark confinement. 
In fact, even if we choose the maximal torus, the area law can be derived. This is 
because the coset G/H is contained in G/H, i.e., G/H C G/H, so that the Wilson 
loop does not feel the whole of G/H, but only feels the components of G/H that 
are contained in G/H. In other words, the variables belonging to G/H — G/H are 
irrelevant for the expectation value of the Wilson loop, as can be seen from the non- 
Abelian Stokes theorem (NAST) that was presented in Ref.[Tl|) and is derived in this 
article. Therefore, a single kind of magnetic monopole is sufficient for confining a 
fundamental quark, since 

7r2(5C/(A)/[/(A - 1)) = 7ri(?7(l)) = Z. (1.3) 

Our results show that two partial gauge fixings 5C/(3) U{2) and SU{3) U{1) x 
C/(l) lead to the same result for confinement as far as the fundamental quarks are 
concerned .O 

The NAST plays the crucial role in this article. The NAST has a number of 
versions which have been derived by many authors. The version of NAST derived 
in this article is based on the idea of Dyakonov and Petrov,El'* who derived an SU{2) 
version and suggested a method of generalization. We derivea version of NAST based 
on the coherent state representation c3) on the flag space, S'lHa* not on the method 
suggested by them. The coherent state representation is used in a different fashion to 

* The possibility of a discrete subgroup has been extensively investigated recently from the 
viewpoint of the Abelian gauge, e.g., the center Zn for SU{N) (see e.g. Ref.|l|)). 
** See Ref.lLQl) for a result of the simulation on a lattice. 
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derive an SU (2) version of the NAST in Refjl|), but the extension to SU{N), N>3, 
was a non-trivial issue which prevented us from presenting immediate pubhcation 
of general SU {N) results. The NAST is not only mathematically (or technicahy) 
important but also physically interesting as we now discuss. 

First, the NAST enables us to write the non-Abelian Wilson loop, 



W'^IA] := ^tr 



V exp ^ig <j> A 



:i.4) 



in terms of the Abelian-field strength (curvature two-form) f = da with the Abelian 
gauge potential (connection one-form) a: 

W^[A] = J [df,{V)]cexp {igf^a^ = J [d^(^)]cexp {^9 Jj^ • (1-5) 

Combining this fact with the_Abelian-projected effective gauge theoriL^(APEGT) 
derived by one of the authors, £3^ we can explain the Abelian dominance E^'E^) in the 
Wilson loop. The APEGT is an Abelian gauge theory obtained by integrating out 
the massive degrees of freedom, i.e., the off-diagonal gluon gauge field with mass 
niji. Hence, the APEGT can be written in terms of the diagonal massless gauge field 
and the anti-symmetric (Abelian) tensor field B^^, together with the ghost and 
anti-ghost fields C and C", where the index i denotes the diagonal components and 
o the off-diagonal ones. Therefore the APEGT is regarded as a low-energy effective 
theory (LEET) which is valid in the long-distance (or low-energy) region R > rrij^ . 
The Abelian gauge field 6^ obtained after the Hodge decomposition of i?^^ can be 
identified with the Abelian gauge field dual to a^. In fact, we can obtain the theory 
with an action S\b\ written in terms of 6^ alone by integrating out all the fields 
other than 6^ in APEGT, and the theory can be rewritten into the dual Ginzburg- 
Landau theory, i.e., the dual Abelian Higgs theory, provided that magnetic monopole 
condensation occurs, i.e., (/c^A:^) / 0. In the dual Ginzburg-Landau theory, the 
coupling constant g in the original theory is replaced by the inverse coupling constant 
1/g which is proportional to the magnetic charge. Therefore, the dual theory can 
be identified with the magnetic theory. On the other hand, the theory with an 
action S[a] written in terms of a* alone is an Abelian gauge theory, but the scale 
dependence of the coupling constant g{ii) is the same as that in the original Yang- 
Mills theory. Thus the low-energy effective Abelian gauge theory exhibits asymptotic 
freedom reproducing the original renormalization-group beta function I3{g). This is 
a manifestation of an approximate weak-strong or electro-magnetic duality between 
two low-energy effective theories described by S[a\ and S\b]. 

Next, the NAST is able to separate the piece oj which corresponds to the mag- 
netic monopole in the Abelian field a = C + u). Indeed, we can write the S'C/(A^) 
version of the 't Hooft-Polyakov tensor describing the magnetic monopole. c3) Hence 
we can separate the contribution of the magnetic monopole in the area law of the 
Wilson loop and explain the magnetic monopole dominance in quark confinement. 
In fact, our derivation of the area law estimates only the monopole contribution, 
Qk = duj. Moreover, the NAST tells us that the essential ingredient in the area 
law lies in the geometric phase, which is concerned with the holonomy group of 
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U {N — 1). Thus we see that quark confinement is intimately related to the geometry 
of Yang- Mills gauge theory, in sharp contrast with the conventional wisdom. 

We present two methods to derive the area law of the Wilson loop by making 
use of the NAST. The first method is to use the APEGT for estimating the diagonal 
Wilson loop; for a sufficiently large Wilson loop (|C| ^ itt-a^), the expectation value 
of the non-Abelian Wilson loop in Yang-Mills theory is reduced to that of the Abelian 
Wilson loop in APEGT: 

W^[A]) = ( exp (ig<f a)) ^ / exp (ig <f a)) . (1.6) 

/ YM \ \ JC J / YM \ \ JC J / APEGT 



Then we can apply the result of Ref.l4), confinement in the Abelian gauge theory, 
to show quark confinement in Yang- Mills theory. 

The second method is to treat the non-Abelian gauge theory directly, without 
going through the effective Abelian gauge theory, based on the novel reformulation of 
Yang- Mills theory in the MA gauge which was proposed by one of the authors. Ej) The 
novel reformulation regards the Yang-Mills theory as the perturbative deformation 
of a topological quantum field theory (TQFT). An advantage of this reformulation in 
the MA gauge is that the derivation of the area law of the non-Abelian Wilson loop 
in the four-dimensional SU{N) Yang- Mills theory is reduced to that of the diagonal 
(Abelian) Wilson loop in the two-dimensional coset (G/H) non-linear sigma (NLS) 
model, at least when the Wilson loop is planar. Therefore the four-dimensional 
problem is reduced to a two-dimensional problem. This dimensional reduction is a 
remarkable feature of the modified MA gauge □ caused by hidden supersymmetry. 
The Yang-Mills coupling constant g of the four-dimensional Yang-Mills theory is 
mapped into the coupling constant in the two-dimensional NLS model. Hence the 
coupling constant is expected to run in the same way as in the original Yang-Mills 
theory, since the perturbative deformation part provides the necessary running, as 
is well known from the loop calculation. For the fundamental quark, we are allowed 
to restrict the flag space -F/v_i := SU{N)/U{1)^~^ to the complex projective space 
CP^-i := SU{N)/U{N - 1). This greatly simplifies the actual treatment. 

Another advantage of this reduction is that the magnetic monopole contribu- 
tion to the Wilson loop in the four-dimensional Yang-Mills theory in the MA gauge 
is equal to the instanton contribution in the corresponding two-dimensional NLS 
model. Indeed, the diagonal Wilson loop can be written as the area integral of the 
instanton density ilx over the area S bounded by the loop C. This correspondence 
may shed more light on the strong correlation between magnetic monopoles and 
instantons observed in the Monte Carlo simulations, since the two-dimensional in- 
stanton is identified as a subclass of the four-dimensional Yang-Mills instanton (see 
e.g. Ref.|ir 



In this article the expectation value of the Wilson loop is estimated by combining 



the instanton calculus and the large N expansion. (See Refs.^8|),|2^), ^0|),|3^) for 
reviews of the large expansion.) We focus on the CP^~^ model corresponding to 
the fundamental quark. First, the instanton calculus is performed within the dilute 



* We must modify the MA gauge slightly in order to keep the supersymmetry, where the 
supersymmetry is expressed by the orthosymplectic group 05*^(412). 
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gas approximation. It is shown that the calculation in the SU{N) case reduces 
to that in the SU (2) case. It is well known that the large N expansion is a non- 
perturbative technique which can be systematically improved. We derive the area 
law to leading order in the large expansion, namely, in the region of large 
and weak coupling g. We hope that our derivation of quark confinement based on 
the dimensional reduction and the large expansion may shed more light on the 
relationship between QCD and string theory, as first suggested by 't Hooft.E^ 

This article is organized as follows. In the first half, we give a derivation of 
the NAST and discuss its implications. Sections 2 and 3 are preparations for § 
4. In § 2 we review the construct of the coherent state on the flag space for the 
general compact semi-simple group G. In § 3, we present the explicit form of the 
coherent state on the flag space for G = SU{N). We define the maximal stability 
group H, which is very important in the following discussion. In § 4, making use 
of the results of §§ 2 and 3, we derive a new version of the non-Abelian Stokes 
theorem for G = SU{N). Although we discuss only the case of SU{N) explicitly, 
it is straightforward to extend this theorem to an arbitrary compact semi-simple 
group G. This version of the non-Abelian Stokes theorem is very interesting not 
only from the mathematical but also from the physical point of view, since the non- 
Abelian Wilson loop is expressed as the surface integral of the two- form (i.e., the 
generalized 't Hooft-Polyakov tensor), which leads to the magnetic monopole. This 
fact is intimately related with the Abelian and magnetic monopole dominance in 
quark confinement, as discussed in subsequent sections. 

In the second half, we derive the area law of the Wilson loop. In § 5 we discuss 
the magnetic monopole in SU{N) Yang-Mills theory. In order to specify the type of 
possible magnetic monopoles, it turns out that the maximal stability group is more 
important than the maximal torus group H. In § 6 Abelian dominance in the Wilson 
loop is shown in the SU (N) Yang-Mills theory in the maximal Abelian gauge based 
on the Abelian-projected effective gauge theory and the non-Abelian Stokes theorem. 
In § 7 we briefly review a novel reformulation of the Yang-Mills theory which has been 
proposed by one of the authors EJ) to derive quark confinement. This reformulation 
is called the (perturbative) deformation of the topological quantum field theory. We 
apply this reformulation to derive the area law of the Wilson loop in SU (N) Yang- 
Mills theory in §§ 8 and 9. In § 9 we show within this reformulation that the area 
law of the Abelian Wilson loop in the two-dimensional nonlinear sigma model for 
the flag space G/H \s sufficient to derive the area law of the four-dimensional Yang- 
Mills theory in the maximal Abelian gauge. At the same time, this derivation leads 
to the magnetic monopole dominance in the area law. In § 8 we demonstrate the 
area law of the Wilson loop in the nonlinear sigma model in an approach based on 
naive instanton calculus. For the fundamental quark, we have only to deal with the 
CpN-^ model. In § 9 we derive the area law based on the large expansion. These 
results imply the area law of the non-Abelian Wilson loop in the four-dimensional 
SU {N) Yang-Mills theory. The final section contains the conclusion of this article. 

In Appendix A, we give derivations of the inner product of the coherent states 
and the invariant measure on the flag space, which are presented in § 3. In Appendix 
B we explain the method of obtaining GP^ and CP^ by gluing the complex planes. In 
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Appendix C, we explain two ways to characterize the element of the flag space and the 
manner of formulating the NLS model using these parameterizations. In Appendix 
D we summarize the large N expansion of CP^^^. In Appendix E supplementary 
material on the expansion is presented. 

§2. Coherent state and maximal stability group 

First, we construct the coherent state \S,,A) corresponding to the coset represen- 
tatives ^ G G/H . We follow the method of Feng, Gilmore and Zhang. c3) For inputs, 
we prepare the following: 
(a) the gauge groupO G and the Lie algebra ^ of G with the generators {T"^}, which 
obey the commutation relations 

[r^,T^]=i/^^cT^, (2.1) 

where the f^^c E'-re the structure constants of the Lie algebra. If the Lie 
algebra is semi-simple, it is more convenient to rewrite the Lie algebra in terms 
of the Cartan basis {Hi, Ea, E^a}- There are two types of basic operators in 
the Cartan basis. Hi and E^- The operators Hi may be taken as diagonal, while 
Ea are the off-diagonal shift operators. They obey the commutation relations 

[H„Hj] = 0, (2.2) 



[Hi, Ea] — aiEa, 
[Ea, E-a] = Cl Hi, 



(2.3) 
(2.4) 



rp TP ^ _ j ^a^pEa+p-, (a + P £ R) , . 

[^-^^]-io, (« + /?0i?,a + /3/O) ^^-^^ 

where R is the root system, i.e., a set of root vectors {ai, ■ ■ ■ ,ar}, with r the 
rank of G. 

(b) The Hilbert space is a carrier (the representation space) of the unitary 
irreducible representation of G. 

(c) We use a reference state within the Hilbert space V^, which can be nor- 
malized to unity: {A\A) = 1. 

We define the maximal stability subgroup {isotropy subgroup) H as a subgroup 
of G that consists of all the group elements h that leave the reference state 
invariant up to a phase factor: 

h\A) = \A)e'^^''\he H. (2.6) 



* Note that any compact semi-simple Lie group is a direct product of compact simple Lie group. 
Therefore, it is sufficient to consider the case of a compact simple Lie group. In the following we 
assume that G is a compact simple Lie group, i.e., a compact Lie group with no closed connected 
invariant subgroup. 
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The phase factor is unimportant in the fohowing discussion because we consider the 
expectation value of operators in the coherent state. Let H be the Cartan subgroup 
of G, i.e., the maximal commutative semi-simple subgroup in G, and Let 7i be the 
Cartan subalgebra in Q, i.e., the Lie algebra for the group H. The maximal stability 
subgroup H includes the Cartan subgroup H = U{lY , i.e., H C H. 

For every element g £ G, there is a unique decomposition of g into a product of 
two group elements, 

g = ^h, i&G/H, heH, (2.7) 

for g G G. We can obtain a unique coset space G/H for a given \A). The action of 
arbitrary group element g £ G on \ A) is given by 

g\A) = ^h\A) = C\A)e''^^''\ (2.8) 

The coherent state is constructed as \(,, A) = (^\A). This definition of the coherent 
state is in one-to-one correspondence with the coset space G/H and the coherent 
states preserve all the algebraic and topological properties of the coset space G/H. 

If r^(a) is Hermitian, then Hj = Hi and El = E^a- Every group element 
g £ G can be written as the exponential of a complex linear combination of diagonal 
operators Hi and off-diagonal shift operators E^- Let \A) be the highest-weight 
state, i.e., Hj\A) = Aj\A), Ea\A) = for a £ R+, where R+{R^) is a subsystem of 
positive (negative) roots. Then the coherent state is given by 123 

\A), rjp e C, (2.9) 

such that the following hold: 

(i) \A) is annihilated by all the (off-diagonal) shift-up operators Ea with a S 
Ea\A) = 0{a G R+y, 

(ii) \A) is mapped into itself by all diagonal operators Hi, Hi\A) = Ai\A); 

(iii) is annihilated by some shift-down operators Ea with a S R-, not by other 
Ep with /3 G Ea\A) = O(some a G i?_); Ep\A) = \A + /3)(some /3 G i?_); 

and the sum is restricted to those shift operators Ej^ which obey (iii). 
The coherent states are normalized to unity: 

{C,A\C,A) = 1. (2.10) 

The coherent state spans the entire space V^. However, the coherent states are 
non-orthogonal: 

{e,A\^,A)^0. (2.11) 

By making use of the the group- invariant measure (i/x(^) of G which is appropriately 
normalized, we obtain 

I \tA)df,im,^\=I, (2.12) 

which shows that the coherent states are complete, but in fact over-complete. This 
resolution of identity is very important to obtain the path integral formula of the 
Wilson loop given in § 4. 



\^,A)=^\A)=exp 



J2 {rilsEfs - fifsE^ 
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The coherent states 1^^,^) are in one-to-one correspondence with the coset rep- 
resentatives ^ € G/H: 

\i,A)^G/H. (2.13) 
In other words, and ^ G G/H are topologically equivalent. 

§3. Flag space and coherent state for SU{N) 

3.1. SU{2) coherent state 

In the case of SU (2), the maximal stability group agrees with the maximal torus 
group U{1) irrespective of the representation. The SU{2) case is well known (see, 
e.g., Ref.^)). The weight and root diagrams are given in Fig.l. 



^ 1- ^ ^ Hi 

— oc cc 

< U--f--l ^ ► 

— oc 

Fig. 1. Root diagram and weight diagram of the fundamental representation of SU{2) where A is 
the highest weight of the fundamental representation. 

The coherent state for Fi := SU{2)/U{\) is obtained as 

\jM = a^)\j,-j) = e^'^-~^'-\j.-j) = (3-1) 

where \ j, —j) is the lowest state, \j,m = —j), of |j, m), and 

J+ = Ji+iJ2, J- = ji, w= (3.2) 

Id cos Id 

Note that (1 + |ti;p)~-^ is a normalization factor that ensures {j,w\j,w) = 1, which is 
obtained from the Baker-Campbell-Hausdorff (BCH) formulas. The invariant mea- 
sure is given by 



2j + 1 dwdw 
Att (1 + \w\ 



For = ^cr^(^ = 1, 2, 3) with Pauli matrices a^^ we obtain J_(_ = ( || ^ ) , and 



1 W 





""^^ e Fi = CP^ = SU{2)/U{1)^ S"^. (3.4) 
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The complex variable w is a CP^ variable written as w = e*''^cot|, in terms 
of the polar coordinate on 5^ or Euler angles, see Ref.l3). We introduce the 0(3) 
vector n as 

n"*^ := sin0 cos 0, := sin0sin0, n^:=cos0. (3-5) 

The relation 

n^{x) = ^a{x)ai(t>b{x) (a, 6 = 1,2) (3.6) 

is equivalent to 

ni = 23f?(0i,^2), n2 = 29(0102), = |</.ip - |02|'. (3.7) 

The complex coordinate w obtained by the stereographic projection from the north 
pole is identical to the inhomogeneous local coordinates of CP^ when 02 7^ 0, 

W - ^(1) + - ^1 + _ 20102 _ 01 , . 

^ " l-na "(|0iP + |02n-(|0iP-|02n"02- ^^-^^ 
The stereographic projection from the south pole leads to 

w = '^^±^=(^X (3.9) 



l + ns V'Pi 

if 01 7^ 0. The variable w is U{\) gauge invariant. Another representation of n is 
obtained by using the parameterization ( |3.4[) of the Fi variable 



n^ = (^|e(u;)V^C(x.)|^) = (0i 0)(J (^-lO) 



This leads to 



^^1 = 1011 {w + w), n2 = —i\4>i\ (w — w), ^3 = |0i| (l — ww). (3.11) 



Indeed, this agrees with (|3.7D if w = (^)*. The entire space of Fi is covered by two 
charts, 

CP' = UiUU2, C/a = {(01,02) GCP1;0„/O}. (3.12) 

3.2. 5*^7(3) coherent state 

For concreteness, we first focus on the SU{3) case. The general SU (N) case 
will be discussed in the final part of this section. The highest weight A of the 
representation specified by the Dynkin index [m, n] can be written as 

A = mhi + n/i2, (3.13) 

where m and n are non-negative integers for the highest weight and hi and /12 are 
the highest weights of the two fundamental representations of SU (3) corresponding 
to [1,0] and [0,1], respectively (see Fig. 2) 
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Therefore, we obtainu 

The generators of SU(3) in the Cartan basis are written as {Hi, H2, E^, Ejs, Ea+f3, 
E-a, E^I3, E-a-fs}, where a and f3 are the two simple roots. (See Fig. 3 for the 
exphcit choice.) 



H2 H2 




Fig. 2. The weight diagram and root vectors required to define the coherent state in the funda- 
mental representations [1,0] = 3, [0, 1] — 3* of SU{3) where A — Fn = ly^ :— (|, is the 
highest weight of the fundamental representation, and the other weights are := ( — ^, 
and u-" ■- (0,-^). 

If mn = 0, (m = or n = 0), the maximal stability group H is given hj H = 
[7(2) with generators {Hi, H2, Ejj, E^^} (case (I)). Such a degenerate case occurs 
when the highest- weight vector A is orthogonal to some root vectors (see Fig. 2). If 
mn 7^ (m 7^ and n ^ 0), H is the maximal torus group H = U{1) x L'^(l) with 
generators {Hi,H2} (case (II)). This is a non-degenerate case (see Fig.4). Therefore, 
for the highest weight A in case (I), the coset G/H is given by 

SU{3)/U{2) = SU{3)/{SU{2) x f/(l)) = CP^, (3.16) 

whereas in case (II), 

SU{3)/{U{l)xU{l))=F2. (3.17) 

Here, CP" is the complex projective space and Fn is the flag space. Therefore, the 
two fundamental representations belong to case (I), and hence the maximal stability 
group is U (2), rather than the maximal torus group U{1) x U{1). The implications of 
this fact for the mechanism of quark confinement is discussed in subsequent sections. 
The coherent state for F2 = SU{3)/U{1)'^ is given by 

\^,A)=aw)\A) ■.= V\w)\A), (3.18) 



* This choice of ^2 is different from that in Refjlip It is adopted so as to obtain the SU (3) case 
when considering the N — 3 case of SU (N) case studied in the next subsecton. 
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H2 



\ 

\ 

\ 

\ 


4 

1 

1 

1 


^< 1 


\ — 


- a(l) / 


\ a(l) 


/ 


\ 


/ 

>^ 


\ 


- a(2) 


-a(3) 



Hi 



Fig. 3. The root diagram of S[/(3), where positive root vectors are given by q'^' = (1,0), o?'^-' = 
(i, -^), and q:^"^-' = (-5^, -^)- Here we have used the same weight ordering as in the SU{N) case 
(see (3.45)) in defining the simple roots. Then the two simple roots are given by := a = a^^\ 



and a := /3 = a 



(3) 



H2 




Hi 



Fig. 4. The weight diagram and root vectors required to define the coherent state in the adjoint 
representation [1,1] = 8 of SU{3), where A = (|, ■^) is the highest weight of the adjoint 
representation. 



with the highest- (lowest-) weight state i.e., 



\C,A)=exp 



\A), (3.19) 
where e~2^ is the normahzation factor obtained from the Kahler potential (ex- 
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plained below): 

K{w, w) := \n[{Ai{w, Jv))"'{A2{w, w))"], (3.20) 
Ai{w,iD) := 1 + + \w2\'^,A2{w,iu) := 1 + |w3p + \w2 - wiW3p(3.21) 

The coherent state A) is normalized, so that A) = 1. We show in Appendix 

A that the inner product is given by 

iC', A\C, A) = e-f^('"''^')e-^[i^(«''V)+^{«','i*)]^ (3.22) 

where 

K{w,w') := ln[l+iDi'wi+tD2't(;2]™[l+«^3'^i'3 + (it'2'-^i'i''tt)3')(^2-tt'i«'3)]"- (3.23) 

Note that K^w, id') reduces to the Kahler potential K{w, w) when w' = w, in agree- 
ment with the normalization. It follows from the general formula (see discussion of 
the SU{N) case) that the SU{3) invariant measure is given (up to a constant factor) 

by 

3 

dfiiO = dfi{w,tv) = D{m,n)[{Ai)"'{A2Tr^ Yl dw^dw^, (3.24) 

a=l 

where D{m, n) = ^(m + l)(n + l)(m + n + 2) is the dimension of the representation. 
For the choice of shift-up (E^i) or shift-down (E^i) operators 

^^'-^TT' ^^'-^VT' ^^'-^VT^ ^'-''^ 

with the Gell-Mann matrices Xa{A = 1, ■ ■ ■ ,8), we obtain 





" 3 




Wl 


W2\ 


exp 






1 






.i=l 
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G F2 = SU{3)/U{lf, (3.26) 



where we have used the abbreviation E±i = E_^^(i) {i = 1,2,3). These two sets of 
three complex variables are related as (see Appendix A) 

n T2 , nT3 T3 ,„ „„x 

-1 = ^, -2 = ^ + ^, -3 = ^, (3.27) 



V2wi, T2 = V2(w2- ) , T3 = V2W3. (3.28) 



or conversely 

n 

The complex projective space CP^ is covered by three complex planes C through 
holomorphic mapsB* (see Appendix B). The parameterization of SU{3) in terms of 
eight angles is also possible in SU{3), just as SU{2) is parameterized by three Euler 
angles (see Ref.f 
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3.3. SUiN) case 

For SU{N) = SU{n + 1), the flag spaceli is defined by 

Fn = SU{n + l)/U{l)'' 3V. 



(3.29) 



We use V to denote an element of -F„ in this definition. Fn is a compact Kahler 
manifold, H*'® which is a homogeneous but nonsymmetric manifold of dimension 
dimc-Fn = ri{n + l)/2. 

Since the fiag manifold Fn is a Kahler manifold, Ei^'H it possesses complex local 
coordinates Wa, a Hermitian Riemannian metric, 



and a corresponding two-form, called the Kahler formJ*l 
which is closed, i.e., 

dQx = 0. 



(3.30) 



(3.31) 



(3.32) 



Any closed form ilx is locally exact {Hk = duj), due to Poincare's lemma. The 
condition (|3.32|) is equivalent to 



7/3 



or 



(3.33) 



This holds if and only if the metric gai3 can be obtained from a real scalar function 
K as 

d d 



(3.34) 



where K = K{w, w) is called the Kahler potential. Then the Kahler two-form is 



obtained from ( |3.31 



as 



K 



iddK. 



(3.35) 



On the flag space, there transitively act two groups, G = SU{n + 1) and its 
complexification = SL{n + 1, C). Any element of Fn can be written as an upper 
triangular (n + l)x(n+l) matrix, whose main diagonal elements are all 1 and whose 
upper n(n + l)/2 elements are complex numbers, Wa £ C: 



f 1 Wi W2 

1 Wn+l ■ ■ ■ 
1 W2n 





\0 



Wn \ 

WSn-S 



1 Wn{n+l)/2 

1 



(3.36) 



* The imaginary unit i is needed to make the Kahler two- form real, since g^fj = ga/i ~ g/Sa- 
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Therefore, we can write 

Fn = SL{n + 1, C)/S_ 9 e, (3.37) 

where B^{B^) is the Borel subgroup, i.e., the group of lower (upper) triangular 
matrices with determinant equal to 1 (Iwasawa decomposition). This definition 
( |3.37] ) should be compared with the first definition ( |3.29| ). The mapping G/H 
C^/B- is a generalization of the stereographic projection in the G = SU{2) ca.se.t3 
The action of the group SL{n + 1, C) on Fn, g : V ^ Vg can be found through the 
Gauss decomposition, 

V-g = TDVg, g £ SL{n + l,C), T£Z_{n + l), Vg£Z+{n + l), (3.38) 

where Z^{n + l)(Z_(n + 1)) is the set of upper (lower) triangular matrices whose 
main diagonal elements are all 1 and D is a diagonal matrix with determinant equal 
to 1 J?he elements of the factors T, D and Vg are rational functions of the elements 
of 5.0 

The group G = SU{N) has rank A^ — 1, and the Cartan subalgebra is constructed 
from (A^ — 1) diagonal generators Hi. Hence, there are A^(A^ — 1) off-diagonal shift 
operators E^, since dim5;7(iV) := iV^ - 1 = (^N -1) + N{N -1). Therefore, the total 
number of roots is N{N — 1), of which there are — 1 simple roots. Other roots 
are constructed as linear combination of the simple roots. Also, there are A^ weight 
vectors. An element of SU{N) is represented by the N x N unitary matrices with 
determinant 1 that are generated by traceless Hermitian matrices, A^^ — 1 linearly 
independent generators T^{A = 1, • • • , A^^ — 1). The generators are normalized as 

tr(r^r^) = ^6ab. (3.40) 

Each off-diagonal generator Ea has a single non-zero element l/\/2. The diagonal 
generator Hm is defined by 

m 

{Hm)ab = , =^W . ^akhk — 'm^a,m+lh,m+l) (3.41) 

^2m(m + 1) 

^ diag(l,---,l,-m,0,---,0). (3.42) 



\j2m{m + 1 



For m = 1 to A'^— 1, the first m diagonal elements (beginning from the upper left-hand 
corner) of H^. are 1, the next one is — m, and the rest of the diagonal elements are 
0. Thus Hm is traceless. The weight vectors (eigenvectors of all Hi: Hj\v) = 
of the fundamental representation N (A^-dimensional irreducible representation of 



* Note that ^ is not necessarily unitary as a matrix under this definition. 
** For n = 1, 



Hence w is the complex one-dimensional representation of 5L(2, C). 
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SU{N)) are given 



11 1 1 



^2' 2^/3' ' V2m(m + 1)' ' ^W^^W / ' 
' 1 1 1 1 \ 

11 1 \ 



\/3'2V6' ' ^2{N - 1)N 



1.™+! = (0,0,---,0, 



^2m(m + l)' ' y/2{N - l)N r 



z.^= (0,0,---,0,^=£t^) . (3.43) 

\ ^2{N - l)N ) ^ ' 

All the weight vectors have the same length, and the angles between different weights 
are the same: 

j,i.yi = E^ yi.yi = — L (fori/i). (3.44) 
2N 2N ^ ' ^ ' 

The weights constitute a polygon in the — 1 dimensional space. This implies that 

any weight can be used as the highest weight. A weight will be called positive if its 

last non-zero component is positive. With this definition, the weights satisfy 

> u"^ > ■■■ > . (3.45) 

The simple roots are given by 

a» = z/i-z/^+i (i = 1,...,A^- 1) (3.46) 

Explicitly, we have 

= (l,0,---,0), 



1^-2,^,0,. ..,oj, 

a^=(o,-^,/f,0,...,o). 



1 ^/Z 



•-^2^- VW^- 
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As can be shown from ( |3.44| ), all these roots have length 1, the angles between 
successive roots are the same, 27r/3, and other roots are orthogonal: 

a^-a^ = l, a'-a^ = -^, {j = i±l) (3.48) 
a' ■ = 0. (j ^i,i±l) (3.49) 

This fact is usually expressed by a Dynkin diagram (see Figj5|). 

0-O--0-0 

ai ai aN-2 aNi 

Fig. 5. The Dynkin diagram of SU{N). 

If we choose as the highest-weight A of the fundamental representation N, 
some of the roots are orthogonal to v^. From the above construction, it is easy to see 
that only one simple root ai is not-orthogonal to u^, and that all the other simple 
roots are orthogonal: 

z^^ • a V 0, i/i • = • = • • • = • a^-^ = 0. (3.50) 

Therefore, all the linear combinations constructed from a^, - ■ ■ , a^~^ are also orthog- 
onal to u^. Non-orthogonal roots are obtained only when is included in the linear 
combinations. It is not difficult to show that the total number of non-orthogonal 
roots is 2(iV - 1) , and hence there are N{N -1)-2{N -1) = {N -2){N -1) orthog- 
onal roots. The (A^ — 2)(A^ — 1) shift operators Eq, corresponding to these orthogonal 
roots together with the — 1 Cartan subalgebra Hi constitute the maximal stability 
subgroup H = U{N-1), since {N -2){N -1) + {N -1) = {N-lf = dimU{N-l). 
Thus, for the fundamental representation, the stability subgroup H of SU{N) is 
given hy H = U{N — 1). In order to describe the coset space G/H, we need only 
(A^ — 1) complex numbers, since 

G/H = SU{N)/U{N -1), (3.51) 

and dimG/H = 2{N - 1). We conclude that GiLH=nP^'^, where CP^-^ is the 
(A^ — l)-dimensional complex projective space, e3'Eb''E3*'E3) which is a submanifold 
of the flag manifold -F/v_i. 

The complex projective space CP" is the compact Kahler symmetric space with 
dim^jCP" = n. The SU{n + 1) group can act transitively on this manifold, and this 
manifold can be considered as a factor space: 



CP" = SU{n + l)/{SU{n) x ^7(1)). 



(3.52) 
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An element of CP" can be expressed using the n complex variables wi, - ■ ■ , Wn as 



(I Wi W2 

1 
1 




\0 








1 
••• 1/ 



G CP". 



(3.53) 



The Kahler potential of P„ (n = N — 1) is obtained as follows. Let H be the 
Hermitian matrix defined hy H = VV^ , V ^ Fn- Consider the Gauss decomposition 



H = VV'^ = TDV, 



(3.54) 



where T G Z-{n + 1) and D is a diagonal matrix with determinant equal to 1. H 
has the upper principal minors A£{VV^){£ = 1, • • • ,n + 1), which are equal to the 
following products of the elements 5i of the diagonal matrix D = diag((5i, • • • , (5„): 



Z\i = <5i, A2 = SiS2, 



n+1 



1. 



(3.55) 



Let (^ = 1, • • • , n) be the Dynkin index of SU{n + 1). Then the Kahler potential 
of Fn is given in the form, 



K{w^ w) = ^ diKi{w, w) = ^ di In Ag{w, w) = In 



1=1 



(3.56) 



The function K(w,w') can also be obtained from the Gauss decomposition of VV'^: 



K{w, w') = In 



1=1 



(3.57) 



The SU{n + 1) invariant measure on P„ is written, up to a multiplicative factor, 

asm 



n(n+l)/2 

dfi{V,V) = p{V,V) n dwo^dwo,, 

a=l 



PiV,V) 



(3.58) 
(3.59) 



where de = 1 for all £. The density p of the invariant measure is calculated from 



P = det(5„^) = det 



(3.60) 



19 



The Kahler potential of the CP^ ^ manifold is given by 

N-l 



K = mln + ) = "^^^ + lll^lll^) ' (3-61) 

where 



N-l 

Iw'lll' := E l^*^!'- (3-62) 

a=l 



Hence, the metric reads 

_ _ {l + \\\w\\\'')5^p-w^W(i 

9af3-^ (1 + |||w|||2)2 • ^'^■^'^^ 

The above construction of the coherent state can be extended to an arbitrary 
compact semi-simple Lie group (see Chapter 11 of Perelomovc2P). 

§4. Non-Abelian Stokes theorem 

In this section we derive a new version of the non-Abelian Stokes theorem 
(NAST) based on the coherent state obtained in the previous section. An advantage 
of this version of NAST is that it is possible to separate the magnetic monopole 
contribution in the Wilson loop and that it is very helpful to understand the dual 
superconductor picture of quark confinement in QCD. 

4.1. Non-Abelian Stokes theorem for SU{N) 

We consider the infinitesimal deviation = ^ + (which is sufficient to derive 
the path integral formula). From (3.19), 



= ^\A) = e-^-^('^''^)exp 

we find 



\A), (4.1) 



+ dC, A\^, A) = exp[ia; + 0{{dw)% (4.2) 

io{x) := {A\iC\x)d({x)\A), (4.3) 

where d denotes the exterior derivative 

d := dxf"^ := dx^da. (4.4) 

Then lo is the one-form 

u; = dx''ojf„ ujf, = {A\i^\x)df,^{x)\A). (4.5) 

Here the x dependence of ^ comes from that of w{x) (the local field variable w{x)), 
i.e., ^(x) = ^{w{x)). 
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The exterior derivative is regarded as the operator 



d = d + d = dwa — h dwj3 



dwr 



(4.6) 



where the operators d and d are called the Dolbeault operators. From the inner 
product (3.22), 



K(w,id') - i [K{w',iD')+K{w,w)] 

^ 1 



{^',A\^,A) = e 



we obtain another expression for u) using the Kahler potential K: 



u; = -{d - d)K = - 



dWn 



d 



dWa 



dw 



d 



K. 



(4.7) 



(4.8) 



The Wilson loop operator is defined as the trace of the path-ordered 

exponent along the closed loop C as 



1 



W"^ \A] := -7tr 
TV 



V exp yig j) A 



(4.9) 



where N is the dimension of the representation {J\f = dim(l/j) = tr(l/j)), and A is 
the (Lie-algebra valued) connection one-form 



A{x) = A^{x)T'^dx^' = A^ix)T^. 



(4.10) 



Consider a curve starting from xq and ending at x. We parameterize this curve by 
the parameter t and define 



Wab{t,to) :-- 
where 



Vex.piig [ dxM^(x)) = Vexpfig [ dtA{t)] 



A{t) := Afj,{x)dx^/dt. 
Then the wavefunction defined by 

Mt) = Wab{t,to)Mto) 

is a solution of the Schrodinger equation 



ab 

(4.11) 
(4.12) 

(4.13) 
(4.14) 



ab 



Note that the Wilson loop operator is obtained by taking the trace of ( 4.11 ) for the 
closed loop, say C: 



dimofrep 

W^[A] = -^tr{W{t,0)) = -^ ^ Waait^O). 



(4.15) 



a=l 
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This implies© that it is possible to write the path integral representation of the 
Wilson loop operator if we identify A{t) with the Hamiltonian as 



H(t) = -gA{t) = -gA^{x)dx''/dt. 



(4.16) 



First, we define the path-ordered exponent by discretizing the time interval t into 
N infinitesimal steps and subsequently taking the limit — > cxd, e ^ with A^e = t 
fixed as 



tr <V exp 



ig [ dtA{t) 
Jo 



N-l 



lim tr <^ P TT [1 + iegAiU^ 

k n=0 



(4.17) 



where t„ := ne, e := t/N. For simplicity, we set to = and = t. Next, we use the 

e repla( 

(4.18) 



coherent state following Ref.l5). On the right-hand-side of (4.17), we replace 

the trace with 



^tr(...) = j d^i{iN){iN,A\{---)\iN,A), 
and insert the complete set (resolution of unity) 

1 = J\^n,A)dfi{i^){^n,A\. {n = l,2,---,N-l) 



(4.19) 



Then we obtain 
1 



tr < 7^ exp 



lim 



ig / dtA{t) 
Jo 

dl^{^N){CN, A\[l + iegA{tN-iMN-i,A)dfi{CN~i) 

X {^N~uA\ [1 + iegA{tN~2MN-2, A)dfxi^N-2) 
■■■dfi{^i){^i,A\[l + iegA{tomN,A) 



AT- 



Af- 



lim f[ f df,{U) l[{Cn+i,A\[l + iegA{tnmn,A) 

hm n / MCn) n [l+iegMtn)] l[{Cn+l,A\Cn,A) 

N 

^™ n n / ^M(^n) exp 



N-l 



n=0 



n=0 
N-l 

l[{^n+l,A\^n,A), 

n=0 



(4.20) 



where we have used = (,N and have defined 

{c.n+1, A\U, A) 

Up to O(e^), we find 

Aitn) := (Cn, A\AitnMn, A) + 0{e^) = {A\^itn)Uitnntn)\A) + 0{e^), (4.22) 
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and 



(C„+l,yl|en,^) = {atn),mtn),A)+e{iitn),mtn),A)+0{e^) 

= exp[e{iitn),mtn),A) + 0{e^)] 

= exp[-ie{A\ii{tJatn)\A) + 0(e')] 

= eMHAmtn)^i{tn)\A) + 0{e% (4.23) 

where we have used yl) = 1. Therefore we arrive at the expression 

TV 

W^[A] = hm n / Matn)) 



n=l ' 
N-l 



exp ige J] + i9-^e(tn)^e(tn)]|vl) .(4.24) 

I n=0 J 



Thus we obtain the path integral representation of the Wilson loop, 



VAV^ + -VdV^ 
9 



\A) 



(4.25) 



where [d^{^)\c is the product measure of d^{w{x),w{x)) along the loop C: 

N 



mi)]c:= ^ lim n I M^n). 

n=l 



(4.26) 



Using the (usual) Stokes theorem, §c=ds ^ 
Stokes theorem (NAST): 



Jgduj, we obtain the non-Abelian 



W^[A] = J[df,{0]cexp [igj^ 



1 

— L 

9 



^^(Olcexp ig 



S:dS=C 



1 
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Here we have defined 



and 



n^(x) := {A\^\x)T^^{x)\A), 
uj{x) := {A\i^^x)dC{x)\A), 

Qk '■= duo. 



(4.27) 



(4.28) 
(4.29) 

(4.30) 



Taking into account ([4.8|), we find that this Qk is identical to the Kahler two- form, 
in agreement with the general statement ( 3.35| ), i.e., 



= duo = {d + d)-{d - d)K = iddK, 



(4.31) 



since the identity d^ = leads to 5^ = = d^,dd + dd = 0. Therefore, the second 
term, co or Hk, in the exponent of the NAST is entirely determined from the Kahler 
potential of the flag manifold. 
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For SU{N), the topological part, 



j:= f u;= Hk, (4.32) 

corresponding to the residual U{N — 1) invariance is interpreted as the geometric 
phase of the Wilczek-Zee holonomy, clP just as in the SU{2) case it is interpreted as 

the Berry- Aharonov-Anandan phase for the residual U (1) invariance. The details of 
this point will be given in a subsequent article. c3) 

4.2. Fundamental representation of SU{N) and CP^~^ variable 
Defining 

uja{x) := {vHx)\A))a, {a = l,---,N) (4.33) 

we can write 

n^{x) := {A\V{x)T^V\x)\A) = Coa{x){T%tMx). (4.34) 

In the CP^~^ case, in particular, the highest-weight state is given by a column 
vector, 

/1\ 



1^) 







(4.35) 



and then we can write 

n\x) := {A\Uix)T^U\x)\A) = Mx)iTXbMx), (4-36) 
where U £ SU{N) and 

Mx) := {U\x)\A))a = Uia{x). (4.37) 
Then n'^(x) can be rewritten as 

n^(x) = {U{x)T^U\x))ii. (4.38) 
Note that the CP^~^ variables (j)a{a = !,•••, N) are subject to the constraint 

N 

J2\Mxr = l- (4.39) 

a=l 

This is clearly satisfied by the unitarity of U, J2a=i Uia{x)Uia{x) = {U (x)C/^(x))ii 
1. 

Now, we examine another expression (the adjoint orbit representation), 

n^{x) = 2ti{U\x)nU{x)T^), (4.40) 

or equivalently, 

n(x) := n^(x)T^ = U\x)nU{x). (4.41) 
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Here 7i is defined by 

N'—^ 1 / AT — 1 1 1 

n = A.iH\..., H^-') = Y: A^H^ = 2diag [-j^, "^y ' " " " ' "iV ) ' ^^'^^^ 



i=l 

where we have used (3.42) and 



A' = , ^ (4.43) 
^Jliii + 1) 



For SJ7(3), when the Dynkin index [m,n] = [1,0] or [0,1], ( 4.42 ) reduces to 
r 1,- /2 1 1\ 1 . /II 2\ 

Note that two elements agree with each other. Hence the adjoint orbit cannot cover 
all the flag space F2. This is the CP^ case. It is easy to see that the two definitions 
(4.36) and ( [4.40| ) are equivalent, 

n^{x) = {V(x)T^V\x))xx = 2tT{nU{x)T^U^x)), (4.45) 

since UT'^U^ is traceless.B 

The correspondence between the -F/v-i variables Wa and the CP^~^ variables 
(t)a is given, e.g., by 

4)1 =Wi, (t)2=W2, (t)N-l=WN-l, = 1- (4.46) 

Thus uj is an inhomogeneous coordinate, uJa = 4>a/(t>N = Wa (a = 1, • • • , — 1), by 
definition. In the CP^~^ case (for the fundamental representation), we can perform 
the replacement 

(yl|/(y)|.l) = 2tr[W/([/)] (4.47) 

if f{V) belongs to the Lie algebra of G (and hence f{V) is traceless, i.e., tr/(y) = 0). 
Thus, we obtain another expression for w, 

uj{x) := 2tv[niU{x)dU\x)] = -i2tT[ndU{x)U\x)], (4.48) 
which is a diagonal piece of the Maurer-Cartan one-form, 

i}:=dUU~\ (4.49) 
It turns out that the two-form l7x is the symplectic two-form, 

Q^=du = 2ii{n[U'^dU, U'^dU]) = 2tr(n[dn, dn]). (4.50) 

Our choice of TL is the most economical one (see Ref.^^ for different choices and 
more discussion of the related issues). 



* When [m, n] — [1, 1], on the other hand, Ti. = diag(l, —1, 0), and all the diagonal elements are 
different. Therefore n(x) moves on the whole flag space, F2. 
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From the Kahler potential of CP^ ^ ( ^.61 ) and the relation ( [4.8[ ), the connection 
one-form u) reads 



which is equal to 



1 WadWa - dWaWa ciN 

LO := z , (4-51) 

2 1 + WaWa 



■ WgdWa 

uj:=i—— , 4.52) 

1 + WaWa 



up to the total derivative. By taking the exterior derivative, we obtain 

(1 + \\\w 



nK = d^ = ^^^^4^7^^^^-" A dw^ (4.53) 



which agrees with the metric ( p. 63 ). 

4.3. An implication of the NAST 

The NAST (4.27) implies that the Wilson loop operator is given by 

W^[A] = J mOhexp (igf^a^ = J [dfimcew [^9j^^_^^f^ -(4.54) 

First, a is the connection one-form 

a := n^A^ + -i^ = (vlU^U), (4.55) 
9 

where is obtained as the gauge transformation of ^ by y € -F/v_i: 

:= VAV^ + -YdV'^ = + -^U^. (4.56) 
9 9 

For a quark in the fundamental representation, we can write 

a = 2tv{nA^). (4.57) 

Therefore, the one-form a is equal to the diagonal piece of the gauge-transformed 
potential A^ . This fact is very useful to derive the Abelian dominance in the low- 
energy physics of QCD (see § 6). 

Next, / is the curvature two-form, 

f ■=da = dC + -duj = dC + -QK, (4.58) 
9 9 

where we have defined the one-form 

C := n^A^. (4.59) 

The anti-symmetric tensor f^j^y can be called the generalized 't Hooft-Polyakov tensor 
for the following reasons: (1) it gives a non-vanishing magnetic monopole (current), 
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where only the second term Qk gives a non-trivial contribution; (2) it is invariant 
under the full gauge transformation, although it is an Abelian field strength. These 
facts are demonstrated as follows. 

First, wc characterize the flag space in complex coordinates. More precisely, 
the target space at each space-time point x G is parameterized by the complex 
variables, = w°'{x). The Kahler two- form is rewritten as 

= ig^^d^w'^d^w^dx'' A dx". (4.60) 

On the other hand, 

:= dw = ^(d^uj^ - d^u^)dx^ A dx'' := ^f^dx^" A dx\ (4.61) 
Then the second piece g'^Qx of / can be written as 

= -^{d^uj, - d,u^) = '-g^-pd^w^d,w^ (4.62) 

and hence 

f„v = d^C, - d,C^ + f^,. (4.63) 
The magnetic monopole current is obtained as the divergence of the dual tensor 

kf, ■■= d,*f^,, (4.64) 
where the Hodge dual of f/^i, in four dimensions is defined by 

//ii^ •= 2^ixi'p(jfpcr- (4.65) 

The first piece dC in / does not contribute to the magnetic current, due to the 
Bianchi identity. On the other hand, the second term f}^ in / can lead to a non- 
vanishing magnetic current, as shown shortly. Here it should be remarked that the 
expression for given in terms of the local coordinate Wa leads to a vanishing 
magnetic current. In fact, if the metric is given by the Kahler potential we have 
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— - ^Hvpa 



^ + .^J!^..^.B ] d.w-^d.w-d^w^ = 0, (4.66) 



2g ^ \dw^dw"dwP dufidw^dufi 

where we have used the antisymmetric property of e^ypa under the exchange of v 
and p, and v and a. However, this does not imply the vanishing total magnetic flux 
or magnetic charge. 
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We recall that this situation is similajLto that of the Wu-Yang monopolecHP com- 
pared with the original Dirac monopole.E^' There are two ways to describe the Dirac 
magnetic monopole. One is to use a single vector potential with (line) singularities, 
called the Dirac string, where the singularities are distributed on a semi-infinite line 
going through the origin of the space coordinates. In the absence of singularities, 
the vector potential gives a vanishing magnetic charge, due to the Bianchi identity, 

^ = ^ B • dS = ^ curlA • dS = J ^ divcurlAdF = 0. (4.67) 

Therefore, the singularities must produce a magnetic charge which has the same 
magnitude as the magnetic charge at the origin, but opposite sign. Another way 
is to introduce more than one vector potential to avoid singularities. Each vector 
potential A" is defined in a sub-region Ua of the sphere such that A" is regular in 
each region Ua and the union of the sub-regions covers the whole sphere. Thus, the 
Bianchi identity leads to zero magnetic flux in each sub-region. Note that we can not 
apply the Gauss theorem divcurlA = 0, since Ua is not a closed surface. The total 
magnetic flux is recovered by summing up all the contributions of the differences of 
vector potentials on the boundary -Bq,/? between two regions Ua and Ujs are 

<P = Y,f curlA" ■dS = Y,{ (A" - A^)d£, (4.68) 

where the minus sign follows from the fact that the orientation of the boundary is 
opposite for neighboring regions. The difference is given by the gauge transformation, 
A° - A^ = VAa,p. This recovers the same magnetic flux as in the former case. 

The variable tu" corresponds to A" in the case of the Wu-Yang monopole. There- 
fore, to show the existence of a non-zero magnetic flux, we must specify the method 
of gluing different coordinate patches on the boundary. These subtleties are avoided 
by using a different parameterization. This generalizes the argument given by 't 
Hooft and Polyakov for the SU{2) magnetic monopole. The antisymmetric tensor 
ffj.u given by (4.58) is the SU{N) generalization of the 't Hooft-Polyakov tensor for 
SU{2). In the SU{2) case, a = 2tr(r^^^) for any representation, and the two-form 
/ := da is the Abelian field strength, which is invariant under the SU{2) transfor- 
mation. Hence the two- form / is identically the 't Hooft-Polyakov tensor, 

fA^) := d^{n^{x)A^{x)) - d,{n^{x)Af,{x)) - -rv{x) ■ {d^n{x) x 5,n(x)), (4.69) 

describing the magnetic flux emanating from the magnetic monopole, if we identify 
with the direction of the Higgs field: 

4>^:=cP^M, |</,|:=y^. (4.70) 

The complex coordinate representation reads 

1 1 
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In general, the (curvature) two-form / = d{n'^A'^) + fix in the NAST is the 
Abelian field strength, which is invariant under the full non-Abelian gauge transfor- 
mation of G = SU{N)B 

Uuix) := d^in''{x)A^ix)) - a,(n^(x)^;^(x)) + ^n(x) • [d^n{x),dMx)]- (4-75) 

The invariance of / is obvious from the NAST (4.54), since is gauge invariant 

and the measure [<i/i(.^)]c is also invariant under the G gauge transformation. In the 
case of the fundamental representation, the invariance is easily seen, becaiise it is 
possible to rewrite (4.58) or ( 4.75 ) into the manifestly gauge-invariant formhil 

f^,{x) := 2tr (n{x)J^^,{x) + ^n(x)[Z?^n(x), D,n(x)]) , (4.76) 

where 

J-^^(x) := d^A^{x) - d^A^,{x) - ig[Ai,{x),Au{x)] (4.77) 

and 

Df,n{x) := d^n{x) - ig[Af,{x),n{x)]. (4.78) 
In fact, we obtain the magnetic charge 

d^xko (4.79) 



/ d^x]-eijkdifjk (i,i, A; = 1,2,3) 

J z)3 z g 
J d^ai^eijk^n- [djii,dkn] 
f 1 i 

/ d^x-eab-'ii ■ [dan, dbii] (a, 6 = 1,2) 
2 g 

- I Qk=-- -vrQ, (4.80) 



* The normalization 

tr(r-^r^) = ]^5ab (4.72) 

holds for any group. For SU{2), trlJ^^T'^T'^) = \ieABC. For Sf/(3), tr(r-'^[r^, T^]) = \ifABC, 
while 

tr(r-'^r^r'^) = ^{ifAsc + dABc). (4.73) 

Here we have used T^T^^ = i[r^,T^] + i{r^,r^}, [r^,r^] = ifBCoT" , and {r^,T^} = 
\5abI + cLbcdT^ , where cLabc is completely symmetric in the indices. Furthermore, we find 

tr{T^T^T'^T^) = —SabScD ~ -fABEfcDE + -dABEdcOE + -{fABEdcDE + fcDEdABE)- 
12 o o o 

(4.74) 

** An explicit derivation of this form is given by Hirayama and Ueno in Ref.|^), where the 
corresponding expression in the adjoint representation is also given. 
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where we have used ( [1.50 ) in the last step. Here Q is the integer- valued instanton 
charge in the NLS model in two-dimensional space, 5-2 = r2 u {00} (see § 8). 
The contribution from the magnetic monopole is replaced by the instanton in the 
two-dimensional NLS model. The magnetic charge satisfies the Dirac quantization 
condition, 

gmg = 2t:Q = 2Tm. (n = 0, ^1, ^2, • • •) (4.81) 

4.4. Explicit forms of uo and Qk for SU (3) and SU (2) 
For SU (3) , we find that lo is given by 

Widwi + W2dW2 



u) = im- 



Ai{w,w) 



wsdws + {w2 - wiws){dw2 - widws - w^dwi) 

A2(W, w) 

up to the total derivative. Hence, we obtain 

Ok = dw = im(Z\i)~^[(l -|- \wi\^)dw2 A dw2 — W2Widw2 A divi 
—W2Widwi A dw2 + (1 + \w2\'^)dwi A dwi\ 
+in{A2)~'^[Aidw^ A dw-^ — {wi + wsW2)dw3 A {diD2 — ivsdivi) 
— {Wl + W3'W2){dW2 — w^dwi) A dius 

+{1 + \w3\'^){dw2 - w:idwi){dw2 - w^dwi)]. (4.83) 
The Kahler potential for F2 is given by 

K{w,w)=H{Air{A2T]. (4.84) 

For CP^, it reads 

K{w,w)=H{Airi (4.85) 

which is obtained as a special case of F2 by setting w-i = Q and n = 0. Hence, we 
obtain 

. widwi + W2dw2 ,, 

uj = im — , (4.86) 

A\\w^ w) 

up to the total derivative, and 

f^K = duj = im(Z\i)~^[(l -|- \wi\^)dw2 A dw2 — W2Widw2 A divi 

—W2Widwi A diD2 -|- (1 + \w2\'^)dwi A divi]. (4.87) 

This should be compared with the case Fi = CP^, 

K{w,w) = mln[{l + \wf)], m = 2j. (4.88) 

For SU{2), we reproduce the well-known results, 
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and 

^ • 7 ,_. dw f\dw , ^ 

iiK = igwwdw Adw = im (4.90) 

(1 + |i(;|^)^ 

The explicit form of fix is necessary to carry out the instanton calculus in the 
following. 

§5. Magnetic monopole in SU(N) Yang-Mills theory 

In the dual superconductor picture of quark confinement, the magnetic monopoles 
give the dominant contribution to the area law of the Wilson loop or the string ten- 
sion. Following the 't Hooft argument, the partial gauge fixing G ^ H realizes the 
magnetic monopole in Yang- Mills gauge theory, even in the absence of an elementary 
scalar field. In the conventional approach based on the MA gauge, the residual gauge 
group was chosen to be the maximal torus subgroup H = U (1)^~^ for G = SU{N). 
This choice immediately determines the type of magnetic monopoles. We now re- 
examine this issue. 

We have learned that the magnetic monopole which is responsible for the area 
law of the Wilson loop is determined by the maximal stability group H rather than 
the residual gauge group H. This is a new feature that appears in SU (N) for > 3. 
It seems that this possibility has been overlooked in the lattice community, as far as 
we know. Indeed, this situation occurs only for SU{N) with > 3. Therefore, we 
must distinguish the maximal stability group H from the residual gauge group H. In 
general, the maximal stability group H is larger than the maximal torus subgroup: 
H = C/(l)^~i C H. Hence, the coset space is smaller than in the maximal torus 
case, i.e., G/H C G/H. 

The existence of magnetic monopoles is suggested by the non-trivial Homotopy 
groups i:2{G/H). In case (II), we have 

7r2(F2) = 7r2(5C/(3)/(?7(l) x ^7(1))) = 7ri(C/(l) x U{1)) = Z + Z. (5.1) 

On the other hand, in case (I), i.e., [m,0] or [0,n], we have 

TT2{CP^) = 7T2{SU{3)/U{2)) = 7Ti{U{2)) = 7Ti{SU{2) X [/(I)) = 7Ti{U {!)) = Z. 

(5.2) 

Note that the GP^ ^ model possesses only local U{1) invariance for any N > 2. It is 
this U{1) invariance that corresponds to a single kind of Abelian magnetic monopole 
appearing in case (I). This magnetic monopole may be related to the non- Abelian 
magnetic monopole proposed by Weinberg et al.E^ The explicit solution for the 
magnetic monopole in SU{3) gauge theories are discussed in Ref]4^). 

This situation should be compared with the SU{2) case, where the maximal 
stability group is always given by the maximal torus H = U{1), irrespective of the 
representation. Therefore, the coset is given by 

G/H = SU{2)/U{1) =Fi = GP^ ^ (5.3) 

and 

7r2(5C/(2)/C/(l)) = 7r2(Fi) = ^GP^) = Z, (5.4) 
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for arbitrary representation. 

For SU{N), our results suggest that the fundamental quarks are to be confined 
when the maximal stability group H is given by ^ = U{N — 1) and 

TT2{G/H) = Tr2iSU{N)/UiN - 1)) = TTiiCP^'^) = Z, (5.5) 

while the adjoint quark is related to the maximal torus H = U{1)^~^ and 

7T2{G/H) = TT2{SU{N)/U{lf-') = 7r2{FN-i) = Z^-\ (5.6) 

This observation is in sharp contrast with the conventional treatment, in which the 
{N — 1) species of magnetic monopoles corresponding to the residual maximal torus 
group f/(l)"^~^ of G = SU{N) are taken into account on equal footing. In fact, the 
NAST derived in this article shows that the fundamental quark feels only the U{1) 
that is embedded in the maximal stability group U{N — 1) as a magnetic monopole. 
This is a component along the highest weight. 

§6. Abelian dominance in SU{N) gauge theory 

6.1. APEGT as a low-energy effective theory 

The Abelian dominance in SU (N) Yang-Mills theory can be explained as follows. 

First, we adopt the maximal Abelian (MA) gauge. The MA gauge for SU{N) is 
defined as follows. Consider the Cartan decomposition of A into diagonal (H) and 
off-diagonal (G/H) pieces, 

A{x)=A^{x)T^ = a'{x)H' + A''{x)T'', {A = 1, ■ ■ ■ , - 1). (6.1) 

In particular, for G = SU{3), 

H' = ^, H^ = Y^ r» = y. (a = 1,2,4,5,6,7) (6.2) 
The MA gauge is obtained by minimizing the functional of off-diagonal fields, 

7^ := y d4x^^«(x)^'^«(x) := J d'xtrG/H{A^{x)A''{x)), (6.3) 

under the gauge transformation. Under the infinitesimal gauge transformation A, TZ 
transforms as 

SaTZ = J d^xAi'^dAAl 

= j d^xA'^'^id^A'' + gr^a^Ai + gf-^a'^A^ + gf^^ A^A^) 

= - j d'^xid^A'"' + gr^alAi'^)A'', (6.4) 

since f^^'^ is completely antisymmetric in the indices and Z"*-' = (T* and 
commute) . Therefore, the condition SaTZ = for arbitrary A leads to the differential 
MA gauge given by 

d^A^^^x) - gr'^'aUx)A^^\x) := {D^ia^^^T = 0- (6-5) 



32 



The SU (N) Yang- Mills theory in the MA gaugeS is given by 

^YM = >^YM + OGF+FP, 

Sym = I <fx^J^^,r^^^, 



^GF+FP 



a 



C [DMA^' + -B 



(6.6) 
(6.7) 

(6.8) 



where (5^ is the Becchi-Rouet-Stora-Tyupin (BRST) transformation, 

6bA^{x) = V^[A]C{x) := d^C{x) - ig[A^{x), C(x)], 

5BC{x)=t]^g[C{x),C{x)l 

6bC{x) =iB{x), 
5bB{x) = 0, 

which is nilpotent, i.e., 8^ = 0. The generating functional is given by 
Zym[J] = I [dA\[dC][dC][dB\eMiS'^M +iSj). 



(6.9) 



(6.10) 



Now we proceed to derive the effective Abelian gauge theory in the MA gauge by 
integrating out the off-diagonal gauge fields (together with the ghost and anti-ghost 
fields), ^'^,C°,C"' and B", as done in Refs.|oD,[l|). Then the SU{N) Yang-Mills 
theory can be reduced to the U{\)^~^ Abelian gauge theory, which is written in 
terms of the diagonal fields, a^,C^,C^ and B^, alone: 



Zym[J] = I [da'][dC'][dC'][dB']exp{iS'XPkGT + iSj), 



(6.11) 



exp{iS'/p^^GT + iSj) = I [dA''][dC''][dC''][dB''] exp(i5^°*f + iSj). (6.12) 



where 



In particular, the partition function reads 

Zym [0] = / [da'] [dC'] [d&] [dB'] eMiSTp^GT) ■= Zapegt- (6.13) 



The Abelian gauge theory obtained in this way is called the "Abelian-projected ef- 
fective gauge theory" (APEGT). It has been shown that the APEGT is an Abelian 
gauge theory whose gauge coupling constant g runs according to the same renormalization- 
group beta function as in the original SU (N) Yang- Mills theory. 



'S'arbgtI'^*' C**, C**, -B*] — I d X 



V(^) 



lliV 



5^(/") 5^(/^o) 87r2 /io 



(6.14) 
> 0, (6.15) 



* In order to fix the residual Abelian gauge group H , we add an additional GF + FP term, 
e.g., d^xiSB \& (a^a" + §By . 
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up to the one-loop level.u Hence 



YM 



(6.16) 



[dA] [dC] [dC] [dB] exp(z5^°^f )/[a^] 
= ZxhEGT I [da^][dC^][dC*][rfi?1exp(i5i«|,lt;r)/[«'l (6.17) 



f[a- 



APEGT 



(6.18) 



6.2. Modified MA gauge 

In the MA gauge, it is expectedH* that the ofF-diagonal gauge fields A'^ have 
non-zero mass r?T-^(m^ ^ 0), whereas the diagonal gauge fields remain massless 
{ma = 0). Therefore, the APEGT obtained in this way can be regarded as the 
low-energy effective gauge theory of Yang-Mills theory. In the framework of lattice 
gauge theory, this prediction was confirmed in numerical calculations by Amemiya 
and SuganumaEa) for G = SU{2) and SU{3). In the framework of the continuum 
gauge field theory, on the other hand, it is more efficient to modify the MA gauge 
into the OSp{D\2) invariant form by making use of the BRST 6b and anti-BRST 
6b transformations, as proposed by one of the authors: Ii3 



Sgf+fp[A,C,C,B] := d X iSb^b^^g/h 



^AfjiA^ 



a 



iCC 



(6.19) 



where a corresponds to the gauge fixing parameter. Note that (6.19) is obtained from 
(6.8) by adding the ghost self-interaction terms and by adjusting the parameter for 
the ghost self-interaction term, since 



-A"'A^"' 
2 ^ 



a 



a 



C D^[a]A'' + -B 



jabi ^a^b^i jabc^a^b^c 



(6.20) 



The special case a = —2 is discussed in previous articles. lil^' 111' For G = SU{2), the 
last two terms reduce to liC^C^C^ = —2G^C^C^ in agreement with the previous 
resultlli (See Ref.|l3D for details). 

In the modified MA gauge (6.19), the non-zero mass generation of off-diagonal 
components was demonstrated analytically (at least in the topological sector), using 
dimensional reduction to the two-dimensional coset non-linear (NLS) sigma model 
(see section IV. C of Ref.[l3|)). Quite recently, dynamical mass generation of the off- 
diagonal gluons has been shown to take place due to ghost-anti-ghost condensation 
caused by the attractive quartic ghost interaction (contained in the gauge fixing term 
of the modified MA gauge), which is necessary to maintain renormalizability in four 
dimensions. The off-diagonal gluon mass obtained in this way can be written in terms 



* See Ref.^ for details. The result of Ref.|^[|) for 5(7(2) can be generalized to SU{N) in 
straightforward way, at least at the one-loop levelEd-* . At the two-loop level, this is not trivial. The 
two- loop result will be given in Ref.M). 
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of the intrinsic scale of the gluodynamics, Aqcd, and hence it is a renormahzation- 
group invariant quanity. (For more details, see Ref.^2|).) Therefore, integration 
of the massive off-diagonal gauge fields is interpreted as a step of the Wilsonian 
renormalization group,G and the APEGT can describe low-energy physics on the 
length scale R > . In this sense, the APEGT can be regarded as the low-energy 
effective gauge theory of Yang-Mills theory. 

6.3. Abelian dominance 

With the NAST for SU{N) just derived, Abelian dominance in SU (N) Yang- 
Mills theory is explained as follows, in a manner based on the same argument as 
that in the SU{2) case.0 

The NAST (4.27) implies that the expectation value of the Wilson loop in the 
SU{N) Yang-Mills theory is given by 

{W^[A])ym = J[dfi{V)]c(^exp(^ig j^^a}j^^^, (6.21) 

where the one-form a is written as 

a = dxf'af, = 2tr{nA^) (6.22) 

for a quark in the fundamental representation. Thus, the one-form a is equal to the 
diagonal piece of the gauge-transformed potential , i.e., a component along the 
highest-weight vector of the fundamental representation of SU{N)x^ By applying 
the above result (6.18) to (6.21), we obtain 

W^[A]) = (expfig (f a)\ (6.23) 

/ YM \ \ Jc J / YM 



explig f a) ) (6.24) 
V JC J I APEGT 

= ZaIectJ [da%d&][d&][dB']^xp{iS%%GTW])^^V {^9 j^a^ , (6-25) 

where we have used in the first equality the fact that the the integration over V is 
redundant after having taken the expectation value {■)ym- This implies the Abelian 
dominance for the large Wilson loop in SU (N) gauge theory in the sense that the 



* Rigorously speaking, all high-energy modes should be integrated out in the Wilsonian renor- 
malization group. Hence, we must integrate out the high-energy mode of the diagonal fields, aj,, 
too. The result is the same as the above, at least at the one-loop level. At the two-loop level, we 
must be more careful in dealing with the high-energy mode (see Refjsi])). 

** Abelian dominance in the low-energy region of S'£i(2)-QCD is demonstrated in Ref.^5|) using 
the result in Ref.hd) combined with the SU{2) NAST.EIP'ES' The monopole dominance is derived 
for SU{2) in Ref.]l5|) by showing that the dominant contribution to the area law comes from the 
monopole piece alone, Qk = duj = tr(n[dn, dn]). 

*** Therefore, the component a along the highest- weight is obtained as an appropriate linear 
combination of a'. Other components orthogonal to the highest- weight do not contribute to the 
expectation value of the Wilson loop, since the APEGT (6.14) is an Abelian gauge theory without 
self-interaction among the gauge fields. 
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expectation value of the non-Abelian Wilson loop in Yang-Mills theory is nearly 
equal to (or dominated by) that of the Abelian Wilson loop in the APEGT, where 
the Wilson loop C is so large that the APEGT is valid in that region, i.e., R,T > vn~^ 
for a rectangular Wilson loop with sides R and T. The APEGT is valid in the range 
^QCD < /U ~ < rUA- 

6.4. Monopole dominance and the area law 

In our framework, the Abelian dominance and the monopole dominance are 
understood as implying the first and the second relations below, respectively: 



W^[A]) ^ (expiig i a) ) (6.26) 

YM \ \ Jc J / APEGT 



^{explif u;) ) . (6.27) 

\ \ Jc J / APEGT 

Numerical simulations show that the monopole part obeys the area law and CTAbel 
exhausts the full string tension obtained from the non-Abelian Wilson loop (i.e., 
monopole dominance in the string tension or area law). 



exp U * iv] ) r^exp{-aAbei\S\), (6.28) 

\ Jc J / APEGT 

while (exp {ig §(^a — i §q ^))apegt does not obeys the area law. This result implies 
the area law of the original non-Abelian Wilson loop. 



W^[A]) ~exp(-a|5|) a^GAbei. (6.29) 

YM 



In Ref.l6), the monopole dominance and the area law of the Wilson loop have 
been demonstrated using the APEGT for G = SU{2). Now this scenario can be 
extended to G = SU{N). The important remarks are in order. 

1. The APEGT has a running coupling constant which increases as the relevant 
energy decreases (asymptotic freedom), and thus the APEGT is in the strong 
coupling region in the low-energy regime (see Fig.^). 

2. The Abelian gauge group in APEGT is a compact group embedded in the 
compact non-Abelian gauge group SU{N). It is the compactness that causes 
the phase transition in the APEGT at ac = j, which separates the Coulomb 
(conformal) phase (a < Uc) from the strong coupling phase (a > Oc)- This 
follows from the Berezinski-Kosterlitz-Thouless (BKT) phase transition in the 
two-dimensional 0(2) NLS model obtained through dimensional reduction. 

3. In the low-energy region such that a{fj,) > ac, the APEGT is in the strong 
coupling phase which confines the quark due to vortex condensation. The above 
described strategy is schematically depicted below. 



36 



D-dim. U(l) Gauge Theory 



Covariant Lorentz gauge fixing 



D-dim. Perturbative U(l) 




deform 
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»-dim. TQFT 






Dimensional reduc 


tion 


D-dim. Perturbative U(l) 


(8) 

deform 


(D-2)- 


dim. 0(2) NLSM 






§7. Reformulation of Yang-Mills theory 



In this section we summarize the novel reformulation of the Yang-Mills the- 
ory proposed in Ref.l3) and elaborated in Ref .p!7|) . This material is necessary for 



subsequent sections. 

7.1. Deformation of topological field theory 

We consider the quantization of Yang- Mills theory on the topological background 
field. We decompose the connection A as 

Af^ix) = Q^ix) + Q^ix) (7.1) 

and identify Q with the quantum fluctuation field on the background field Q. For 
arbitrary but fixed background field il, the generating functional is given by 



Z[J, n] := J [dQ][dC][dC][dB] exp{iSYM[^ + Q] 

+iSGF[Q,C,C,B]+i{J^-Q^)}, (7.2) 
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where S'yMlv^] is the usual Yang- Mills action, 

Sym[A] = -J d^'x^J^I^Ml^^'^M, (7.3) 

and Sgf is the gauge-fixing and FP ghost term for the quantum fluctuation field Q: 

Sgf[Q,C,C,B] :=- j d^x ids tTG 

We wish to retain gauge invariance for the background field f2 even after the gauge 
fixing for Q. This is realized by choosing the background field (BGF) gauge fixing 
condition 

F^[Q] := V^^[n]Q^'^ = 0. (7.5) 

In fact, in the BGF gauge, Z[J, f2] is invariant under the gauge transformation of 
the background field; the infinitesimal version is given by 5f2i^ = P^[J7]u;. Hence, the 
theory with the action 

Seff[J, f2] := -ilnZ[J, f2], (7.6) 

is defined only on the space of the gauge orbit. Suppose that the background field 
satisfies the equation ii''^[i7] = 0. In order to consider the quantized Yang- Mills 
theory on all possible background fields satisfying the equations F^[i7] = 0, we 
define the total generating functional 




[dn^][dC][dC][dB]Z[J, n] exp{iSTQFT[^, C, C, B]) exp[i(J^ • ni")] 



= I [df2^][dC][dC][dB] eMiSefflJ, ^] + iSTQFT[fi, C, C, B] + i{J^ ■ ^?'^)^7■7) 

where Stqft[^,C,C,B] corresponds to the gauge-fixing term for the background 
field J7^. In order to describe the magnetic monopole as a topological background 
field in Yang-Mills theory, we choose the MA gauge for f}^, 

STQFT[f2,C,C,B] :=- j d^x i6B tla/H 

where 

:= D^''^[0']0^''' := {8^6"^ - g^^^Q^')Q^'^ {i = I, ■ ■ ■ , N - I). (7.9) 

Note that the trace is taken on the coset G/H, not on the entire G, in the MA gauge. 
Under the identification 

n^{x) := -U{x)d^U\x), Q^(x) := U {x)V ^,{x)U\x) , (7.10) 

we assume that all the topologically non-trivial configurations come from whereas 
V denotes topologically trivial configurations. Therefore, V changes under the small 



C[F[Q] + -B 



(7.4) 



a 



C{F[n] + -B 



(7.8) 



38 



gauge transformation, while Q includes the effect of large gauge transformations. 
Therefore, we must take into account the finite gauge rotation [/, without restriction 
to the infinitesimal gauge transformation. Under the identification (7.10), the Yang- 
Mills action is invariant, 

Sym[A] = Sym[0 +Q]= Sym[V] = J d''x^J^^,[V]J'^'''^[V], (7.11) 

while the gauge-fixing term (7.4) is chang into 



(7.12) 



This implies that the background gauge for Q is changed into the Lorentz gauge for 
V, d^V^ = 0. Then the generating functional in the background field Q is cast into 

Z[J, Q]= j [dV][d-i][d^][d0\ exp [iSYKim + i^GF[V,7,7,/?] + i(^/. • UV^U^)] , 

(7.13) 

where V;i,7,7 and /? are defined by the adjoint rotations of Q^, C, C and B respec- 
tively: 

:= U^Qf,U, 7 := U^CU, 7 := U^du, /3 := U^BU. (7.14) 
Thus the total generating functional reads 

Z[J] = f [dU][dC][dC][dB]Z[J, n] exj>{iSGF[^, C, C, S]) exp[i(J^ • 



[dU][dC][dC][dB] exp{iSeff[J, ^2] + iSrQFTif^, C, C, B] + i{J^ ■ f?^)}, 

(7.15) 

where we have made a change of variable from f2 to U. The measure [dU] is invariant 
under the multiplication 

U{x) U{x)U{x), (7.16) 
which leads to the finite gauge transformation of the background field, 

Q{x) U{x)Q{x)U\x) + -U{x)dW{x). (7.17) 

9 

The modified MA gauge proposed in Ref. El* leads to the TQFT with the action 



StqftI^, C, C, B] := I d^x i^s^Btrc/H 



-Qr.n^' - -iCC 

2 ^ 2 



. (7.18) 



The expectation value of the functional f{A) of A is calculated as follows. If 
the functional is of the form f{A) = g{Vfj,,U)h{U), then the expectation value is 
calculated according too 

{f{A))YM = MV^, U))^YMhiU))^QpT, (7.19) 



* See Remark 10.1. 
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where the sector of the perturbative deformation is defined by 

{{■■■))^YM = Z;^^,J[dV][dj][dj]mexp{tSYM[V]+iSGF[V,j,j,f3]} (•••), 

ZpYM ■■= J [dV] [d-f] m m exp {iSvM [V] + iScF [V, 7, 7, /3] } , (7-20) 
and the sector topological field theory is defined by 

{{■■■))tqft ■■= ZtqftJ [dU][dC][dC][dB]eMiSTQFT[n,C,C,B]}i...), 
Ztqft ■■= J [dU][dC][dC][dB] e^iSTQFTi^, C, C, B]}. (7.21) 

This reformulation of the Yang-Mills theory is called the "perturbative deformation 
of a topological quantum field theory". The expectation value ((•••))pyA/ 
field V is calculated using a perturbation theory in terms of the coupling constant 
g. On the other hand, the expectation value ((• • •))tqft should be calculated in a 
non-perturbative way to incorporate the topological contribution. Here C/ is a com- 
pact gauge variable corresponding to a finite gauge transformation. In the instanton 
calculus, the integration measure [dU] is replaced with the (finite-dimensional) inte- 
gration over the collective coordinates of the instanton. 

7.2. Dimensional reduction to the NLS model in the MA gauge 

It has been shownEl) that, due to the hidden supersymmetry OSp{D\2), the 
TQFT part (7.18) is reduced to the {D — 2)-dimensional coset (G/H) nonlinear 
sigma (NLS) model, 



Snlsm[U,C,C] := an / trc/H 



2 ^ ^ 2 



(7.22) 



where, for the matrix element ilab (see Appendix C) 

"1 
2 



tr, 



G/H 



-n^{x)f2^{x) 



Wx))abi^^^.ix))ab- (7.23) 

a,b:a<b 



By making use of the complex coordinates in the fiag space G/H, the action can be 
rewritten as (see Appendix C) 

where we have omitted to write the decoupled ghost term, C{x)C{x). In particular, 
for D = 4, 

Snlsm = —,j^ dzd-z g^-, + j ' (^-25) 

where z = x + iy = xi + ix2 € C = and dxdy = dxidx2 = ^dzdz. The G = SU{2) 
case is analyzed in Ref.l3), and in that case we have 

an f 1 f dw dw dw dw\ 

Snlsm = ^ dzdz (^^^^ (, 9? + 9? S^J ' ^^-'^^ 
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The above described strategy is schematically depicted below.l 
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§8. Area law of the Wilson loop (I) 



In this section we derive the area law of the Wilson loop based on the instanton 
calculus. 53), 54), 55) A more systematic estimation is given in the next section based 
on the large expansion. 

The static potential V{R) is evaluated from the rectangular Wilson loop C with 
sides T and R according to 



1 



ViR) :=- lim -ln{W'-[A])YM,- 

I — >oo 1 



The (full) string tension a is defined by 



a :- 



lim 



1 



ln(T^^[^])yM4, 



A{C)^oo A{C) 



(8.1) 



•2) 



where A{C) is the minimal area of the surface spanned by the Wilson loop C . Of 
course, the rectangular loop has minimal area: A[C) = TR. 
Using the NAST for the Wilson loop operator. 



W^[A] = J[dn{C)]cexp (^ig ^ n'^A^ + i £ uj^ 

we can write its expectation value in the Yang-Mills theory asli^'S-* 

{W^[A])ym, 



.3) 



* When we encounter the NLS model obtained through dimensional reduction in the following 
sections, the coupling constant g should be replaced as 



9 ^ -g 

a 



(7.27) 



since we do not express the a dependence explicitly. 
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{( 



exp 



ig f dx^'n^{x)V^{x) 



c 



exp 



pYM4 



I (p UJ 



•4) 



TQFTa 



where the measure /[<i/i(.^)]c gives a redundant factor which does not affect the 



string tension in the area law, as shown in Section IV. C of Ref. 15) and hence it is 
omitted in the following. 

8.1. Perturbative expansion and dimensional reduction 

In our reformulation, the field V^{x) is identified with the perturbative defor- 
mation. Thus we expand the first exponential of (8.4) in powers of the coupling 
constant g: 



exp 



ig j dx^'n^{x)V^{x) 



pYMi 



ti Jc ' Jc ' 



dx!f,"n^^{xi)n^'^{x2) ■ • •n^"(2;„) 



C 



Hence, we obtain 

)ym4, 



(8.5) 



= (^exp 



I (f> UJ 



n\ Jc 



c 



X ^n^i (xi)n^2 (-3,2) • • • n^" (x„) exp 

X (V^/ {x,)Vf:^ (X2) • • • V;f: (x„))pyM4 



I (f) UJ 



TQFT4, 



(^exp 



I (f) UJ 



TQFTi 



1 + E 



(i5)" 



n=l 



n! 



c 



dxf i dxf •••idx^(V;/(xi)V^,^(x2)---V;:(x„))pyM4 



(n^i(a;i)re-^2(-3;2) . . . n-^"(xn) exp [i§c^])TQFTi 
(exp [ifc-^D^QFTi 



(8.7) 



We restrict the Wilson loop to a planar loop. This choice has the following advantage. 
For the planar Wilson loop C, the Parisi-Sourlas dimensional reduction leads to the 
following identities: 113* 



(^exp 



I (f) UJ 



= ( exp 

TQFT4 \ 



I (f) UJ 



NLSM2 
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and 



n^'^ (xi) • ■ • n'^"{xn) exp i f to 
\ I Jc . 



n^^{xi) ■ ■ ■ n^"(x„) exp 



TQFTi 



c 



NLSM2 



where xi, ■ ■ ■ , G C C R^. This is because, in the case of the fundamental repre- 
sentation of SU{N), and to can be written in terms of [7 as (see (4.28), (4.29)) 

n^{x) = Uia{x){T\bUib{x) = Mx){T\bUx), (8.10) 



^{x) = ^[Ula{x)dUia{x) - dU-ia{x)Uia{x)] 

% — — 

= -[(t)a{x)d(t)a{x) - d(t)a{x)(t)a{x)]. 



Taking the logarithm of the Wilson loop (8.7) and expanding it in powers of the 
coupling constant, we obtain 

\n{W^[A])YM, (8.12) 



= ln( exp 


i <j) u) 














+ ln 


1 - 


.li 




2 Jc 


= ln( exp 


i <j) u) 









NLSM2 



{n'^{x)n^{y) exp [i Jg Qk\)nlsm2 

(exp [i§ciA)NLSM2 



NLSM2 



2 



(n^(x)ra-^(y)exp [i§c^])NLSM2 



+ 0(/), (8.13) 



(exp [i§c'^\)NLSM2 

where we have defined the two-point function 

Gi^{x,y) := {V^{x)V^{y))pYM,. (8.14) 

In the rest of this section we focus on the first term in (8.13). The remaining 
terms will be estimated in the next section. 

8.2. The instanton in Fat-i and CP^~^ models 

We wish to demonstrate the area law for the expectation value. 



exp \i f oj 



NLSM2 \ 



NLSM2 



= ^NLSM2 J [dK'^, w)] eyip{-SNLSM2 [w, w]) exp f2, 
where the two-dimensional NLS model is defined by the action 

S^r.SM2 = -^J^^dxg^^^^^ 



TT 



g^n) Jc [~d^~df + dz dz 



(8.15) 
(8.16) 

(8.17) 
(8.18) 
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where 



z = X + iy = xi + ix2 G C = R^, dxdy = dxidx2 = ^dzdz, (8.19) 

and g{fi) is the Yang-Mills coupling constantE 
Note that the action satisfies the inequality, 

Snlsm = 2g2(^^^ J^2 '^^^9a0aW°' ± ieabdbW°'){daw'^ ± ieacdcw'^)* 



±^-o7-T 1 d'xeatg^,^^ (8.20) 



IT 



g^ilj) " "'^ dxa dxb 



This inequality is saturated when satisfies the equation daW°' ± ieabdbw"' = 0, 
which is equivalent to the Cauchy-Riemann equation, 

B^w'^ ■= {di + id2)w° = 0. (8.22) 

The solution w"' = /"(-z) is an arbitrary rational function of z. 

The finite action configuration of the coset NLS model is provided by the in- 
stanton solution, which is a solution of the Cauchy-Riemann equation ( |8.22| ). It is 
known© that the integer- valued topological charge Q of the instanton in the -F/v-i 
NLS model is given by the integral of the Kahler 2-form over : 

1 /■ _ r d^x dw°' dvjf^ 



vr J J-R2 TT dXa dXb 



dzdz f dw°' dw^ dw" dw^ \ 

This is a generalization of the Fi case {N = 2), where 

if dwdw i f dzdz f dw dw dw dw\ 

^^2^Jc{l + h|2)2 " 2^ Js2 (1 + \w\^f [~dl~dl ~ 5? ■ ^ ■ ' 

Thus the instanton solution is characterized by the integral topological charge 
Q £ Z. For instanton (Q > 0) and anti-instanton (Q < 0) configurations with a 
topological charge Q, the action has 



Snlsm = ^\Q\. (8.25) 
9^ 



The metric in the Kahler manifold F^-i is obtained in accordance with the 
_ ^ _d_ a 

0/3 dw<^ dwl^' 



relation g^g = jA^-rprpK from the Kahler potential 



N-l N-l 

K{w,w) = ^ deKe{w,w) = ^ de\nAi{w,w), (8.26) 



* The running of the couphng constant is given by the perturbative deformation in four- 
dimensional Yang-Mills theory, as in (6.15). 
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with the Dynkin indices d£{i = 1, - ■ ■ ,N — 1). Then the integral / f^x over the whole 
two-dimensional space reads 



. N-l 

/ = ttQ = vr ^ diQe, (8.27) 



where Qi are integer- valued topological charges. Hence, Qk{x)/t^ is identified with 
the density of the topological charge (up to the weight due to the index d^). 

Now we consider the CP^~^ model. If we identify Wa with the inhomogeneous 
coordinates, e.g., Wa := the metric ( |3.63| ) can be rewritten as 



= ^:n2^2 ^ (8-28) 



(Ml 

where 

N 

ll'AlP:=El'A«l' = l'^^l'(l + IIIHIl')- (8.29) 
«=i 

The action of the CP^~^ model is given by 

ScpN-i = 4 / d^xg^f^{w)d^w''d^wf^, (8.30) 

or equivalently, 

ScpN-i = 4 / d''xg^i^{(P)d^rdf,^. (8.31) 
Under the constraint ||0|P = 1, the action can be written as 



bcpN-i - — 

= ^J d'^xiid^^d^r + ird,r)i^f'd,<p^)]. (8.32) 

This agrees with the action of the CP^~^ model presented in Ref.^3|). (See Appendix 
C for more details). 

8.3. Area law in the dilute instanton-gas approximation 

If the Wilson loop is large compared with the typical size of the instanton, 
in (8.16) counts the number of instantons re^p minus anti-instantons n!!* 
that are contained inside the area 5 C bounded by the loop C: 

f?^ = 7r(n!^ - n!_"). (8.33) 

(See Fig.7.) Thus, the expectation value (exp {i fix)) nlsm is calculated by sum- 
ming over all the possible cases of instanton and anti-instanton configurations (i.e., 
by the integration over the instanton moduli). In this calculation, we use 

2 2 

Snlsm = ^\Q\ = ^(n+ + n„), n± = v}^ + n°^\ (8.34) 
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Fig. 7. Instanton and anti-instanton configuration and the Wilson loop C in a finite-volume region 
Li X L2 in the two-dimensional plane embedded in R'*. 

where n^* (n""*) is the number of instantons (anti-instantons) outside S, and n+ 
(n_) is the total number of instantons (anti-instantons). For the quark in the fun- 
damental representation N (di = 1, d2 = = ■ ■ ■ = d^-i = 0), this is easily carried 
out as follows. 

For the SU{3) case with [1,0], an element ^ G F2 is independent of W3, so that 
W3 is redundant in this case. Hence, it suffices to consider the CP^ model for the 
fundamental quark (up to Weyl symmetry). For CP^, the Kahler two- form is given 
by (4.87) 

f^K = im{Ai)^'^[{l + \wi\'^)dw2 A did2 — W2Widw2 A didi 

—W2Widwi A dw2 + (1 + \w2\'^)dwi A dwi]. (8.35) 

When 102 = 0, Qk reduces to 

=i{l + \wl\^)-'^dwl^diBl. (8.36) 
Similarly, when i^i = 0, we have 

= i{l + \w2\^)~'^dw2^d€J2. (8.37) 
For a polynomial Wa = Wa{z) in z = x + iy order n, we find the instanton charge 

j Qk = t^Q, Q G Z. (8.38) 

This is the same situation as that encountered in SU{2), in which case 

j ^K = 2jTrQ, (8.39) 

where j = 1/2 corresponds to the fundamental representation.©'© Thus the Wil- 
son loop can be estimated by the naive instanton calculus. In fact, the dilute instan- 
ton gas approximation leads to the area law for the Wilson loop (see Ref.|T3|)). Here 
the factor vr is very important. The integral of the Kahler two-form f2x is a multiple 
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of TT. If we had a factor 27r, the area law would not hold, just as in the j = 1 case of 
SU{2). 

For the SU{N) case with Dynkin index, [1,0, • • • ,0], it suffices to consider the 
(jpN-i niodel. When Wa ^ ^ and Wf, = for all h ^ a, the Kahler two-form (4.53) 
for CP^~^ reduces to 

Qk = «(1 + \Wa\^)''^dWa A dWa, (8.40) 

with no summation over o. Therefore the above argument can be applied to SU{N) 
for any N. This implies the confinement of fundamental quarks in SU{N) Yang- 
Mills theory within the approximation of a dilute instanton gas. This naive instanton 
calculation can be improved by including fluctuations from the instanton solutions 
following Ref .p5|) , and this issue will be discussed in detail in a subsequent article. 

§9. Area law of the Wilson loop (II) 



The derivation of the area law of the Wilson loop in the four-dimensional Yang- 
Mills theory in the MA gauge is reduced to demonstrating the area law of the diagonal 
Wilson loop in the two-dimensional coset NLS model. In this section we complete 
a derivation of the area law of the Wilson loop in the fundamental representation. 
Here we use the large N expansionS^' ^)>ii)'il)'i^) for the coset NLS model. (See, 
e.g., Refa.p^),|2^), 30),^) for reviews of the large expansion.) 

To perform the large N expansion, it is convenient to introduce the new variables 



Pabix) := 4>aix)4>bix) = Uia{x)Uibix) 

which are used to rewrite 



n\x)V^{x) = Pabix)Vf{x), 



(9.1) 
(9.2) 

(9.3) 



where ^ = 1, • • • , - l and a, 6 = 1, • • • iV for SU{N). 

In terms of the above variables, the expansion (8.7) of the expectation value of 
the Wilson loop operator in powers of the coupling constant g is rewritten as 



^exp 









+ E 



{igT 



n=l 



c 



TQFTi 



c 



dxf'rr 



x( PaibA'^l)Pa2b2{x2) ■ ■ ■ Pa„b„iXn) exp 



(9.4) 



TQFTa 



^exp 









TQFT4 



n=l 
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(-Pgifei {xi)Pa2b2{x2) ■ ■ ■ Pa„b„(Xn)exp [i§c^])TQFTi 



(9.5) 



The diagrams needed to calculate this expectation value are drawn in Fig. 9 based 
on the Feynmann rule given in Fig. 8. Here it should be remarked that the definition 
of the Wilson loop operator 



1 



{W^[A]) :=-^(tr 



T^exp y9 j> dx'^A^{x) 



(9.6) 



yM4 



includes the normalization factor AA^^ and that the expectation value ( p.(j| ) of the 
Wilson loop may have a well-defined large limit. In particular, in the zero coupling 
limit, the expectation value reduces to one. 



(a) 



(b) 



(c) 



a b 



(d) 



(e) 





Fig. 8. Feynmann rule and the corresponding large A'' rule (double line notation due to 't Hooft) 
in QCD. Propagators: (a) quark propagator, (b) gluon propagator. Vertices; (c) quark-gluon 
vertex {qym), (d) three-gluon vertex (grM), (e) four-gluon vertex [gyM)- 



9.1. Large N expansion and dimensional reduction 

It is known ii) that only the planar diagrams contribute to the expectation value 

i^Z"' i^MZ"' (^2) • • • AlrJ"' iXn))YM, (9.7) 

in the leading order of the large N expansion. See Fig. 9. However, it is extremely 
difficult to sum up the infinite number of terms belonging to the leading order of the 
large expansion and to obtain a closed expression in the four-dimensional case. Of 
course, this does not exclude the possibility that the closed expression obtained by 



48 



summing up all the leading diagrams may exhibit the area law. In fact, this strategy 
has been applied in the two-dimensional case and has successfully lead to the area 
law (see, e.g., Ref.|6ll)). 

For the planar Wilson loop C, we have already shown that the Parisi-Sourlas di- 
mensional reduction occurs and that the TQFT sector reduces to the two-dimensional 
coset NLS model, i.e. the NLS model on the flag space Fn-i. Hence, we obtain 



e Jc 



TQFTa 



NLSM2 



(9.8) 



and 



TQFTa 



I NLSM2 

(9.9) 

where xi,---,Xn G C C R^. For the quark in the fundamental representation of 
SU{N), the relevant NLS model can be restricted to the CP^^^ model. 

We now return to the expression (9.5) obtained by way of the NAST and apply 
the large N expansion to the (perturbative) deformation sector and the TQFT sector 
simultaneously. Taking the logarithm of the Wilson loop, therefore, we obtain 



\n{W'^[A])YM, 
In^ exp 

ln( exp 



+ In 



1 



CP' 



CP' 



f[C] + 0{\^/N'^ 
f[C] + 0{X^/N^), 



where we have defined the two-point correlation function 



(9.10) 



(9.11) 



and 



f[C] ■.= ^f^dx^f^dy'^Gff''{x,y) 



{Pab{x)Pcd{y)ey:p [i§cuj])^p^ 



JV-l 
2 



(exp [i§c^])cp. 



N-l 
2 



(9.12) 



It should be remarked that the f[C] is 0{X/N). This is different from the result 
of the usual large N expansion of the Yang- Mills theory, i.e., dia gr am (a) of Fig. 9, 
which is 0(A). This fact is shown as follows. In a previous article, IIZP it is shown that 
the perturbative sector obeys the Lorentz-type gauge fixing, 9^V^ = 0, by virtue of 
the background gauge. Here we adopt the Feynman gauge to simplify the calculation. 
Then the propagator for reads 



AB i 



47r2 



x-y\ 



(9.13) 
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Then we find 



Gff\x,y) := {Vfix)Vfiy)),YM, 



{V,{x)V^{y))pYMAT'')ab{T^U 

Af2-1 



(9.14) 
(9.15) 

(9.16) 



A=l 



where for G = SU{N), 



iT^)aJ,{T\d = \ (SadSbc - ^<^a6<^cd) ■ (9-17) 



This impUes that 



A=l 



2N 



Sab = C2{R)Sab, 



(9.18) 



where C2{R) is the quadratic (second order) Casimir invariant of the fundamental 
representation. If G = U{N), the relation is simplified as J2A=ii'^^)abiT^)cd = 
^Sadhc- The difference between SU{N) and U{N) disappears in the large N limit. 
To leading order, we can set (see Appendix E) 



{Pab{x)Pcd{y) exp [i §c to]) 



CP' 



= {Pabix)Pcd{y)) 



CP.' 



(9.19) 



Thanks to the SU (N) invariance, it is easy to see that 

{Pabix)Pcdiy))cP^-^ = (^adhc - J^SabScdj Q{x,y), (9.20) 



where 



Q(,x,y) 



This leads to 



{<Pa{x)4>bix)4>biy)4>aiy))cp'^-^ - Tji'Pai^)Mx)My)4'biy))cp'^ 



g'Gff'{x,y) 



{Pab{x)Pcd{y) exp [i §^ Uj])cp, 



(9.21) 



JV-l 

2 



(exp[i/c'^])cp. 



AT-l 
2 



Ar2_i 



= 5^,D{x,y)g^ ^ {T'')ab{T^U i^adhc - ^^ab^d) Q{x,y) 
A=l ^ 

= 5^,„D{x,y)g'^ E 
A=l 



tv{T^T^) - -^tr(r^)tr(r^) 



Q{x,y) 



SnuD{x,y)g 



■'^'-^Q{x,y). 



(9.22) 
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Note that Q{x, x) is 0{N and hence f[C] is 0{X/N). This is because we consider 
the large expansion of the CP^~^ model (see Appendix D) with the Lagrangian 
(see(C-23)) 



jC.(JpN-l 

where V^{x) is the composite gauge field, 



vr 



(9.23) 



under the constraint 



^\x)^{x) := ^''{x)^''{x) 



1. 



(9.24) 



(9.25) 



It is not difficult to show that the above estimation gives the correct order for 
the higher-order terms, e.g., (b) and (e) in Fig. 9, by making use of the relations 



(|4.73| ) and (|4.74|) . For example, 5^(V;f^(xi)V^2(x2)V^^(x„))pyM4 is proportional to 



igifABC^ and iffV^-^^tr[T^r^T^]/iV3 ^ _g4fABCfABC/^^^3-^ ^ o{g^), since 
fABCfABD ^ C2{Adj)6^^ and C2{Adj) = N. 

Another way to understand this result is based on the idea of the reduction 
of degrees of freedom that are responsible to the Wilson loop. The flag space has 
dimension dimF/v-i = A^(A^ — 1), whereas CP^~^ has dimension dimCP^"^ = 
2(A^ — 1). Therefore, the number of relevant degrees of freedom is reduced for 
the fundamental quark for large A^, since dimCP^-i ^ 2dimFAr_i/iV for large N. 
Indeed, this result is expected from the NAST given by (4.54), 



c 



(9.26) 



The Abelian gauge field a = tr(TCA) has only two physical degrees of freedom, while 
the non- Abehan gauge field A = A^T"^ in the Wilson loop (|9j) has 2{N'^ - 1) 
components. Thus the large expansion is reduced to a perturbative expansion in 
the coupling constant g. In this sense, the large A^ expansion combined with the 
NAST justifies the identification of the deformation part with the perturbative part. 

Then we find that the f[C] is of order 0{\/N). Therefore, to leading order in 
the large A^ expansion, the static potential and the string tension are given by 



ViR) 



a 



lim — In/ exp 



lim — — 

R,T^oo RT 




ffl + 0(^' 



A2 



(9.27) 



+ /[C]1 +0(^).(9.28) 



Now we proceed to estimate the second term, f[C]. We will show that the second 
term gives at most the perimeter law, so that the area law (if it exists) is provided 
by the first term. Because of the factor 5^,^ in (V^(x)V,f (y))pyA/4) only integration 
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between parallel sides dx and dy gives a contribution to f[C]. Thus f[C] is reduced 
to 

f[C] = ^ f^dx^' f^dy^D{x,y)Gc{x,y), (9.29) 

where we have defined the correlation function for the composite operators as 

{Pab{x)Pcdiy)ex.];)[i§^Lj])^pN-i 

Gc{x,y) := 2{T^U{TXd 1 ... ' (9-30) 

(exp [i§(j^\)cp^' 



3JV-1 



^ 2{T^U{T^),d [Sad^bc - ^^abScd^j Qix, y) (9.31) 
= {N'-l)Q{x,y). (9.32) 
First, if we restrict our consideration to topologically trivial configurations, i.e., 

n^{x)n^{y) ^ Pab{x)P,d{y){T^)ab{T^)cd = n^(oo)n^(oo) ^ ^(1 - iY"!), 

(9.33) 

then we obtain Gc{x,y) = 1 + 0{N^^) and 

f[C] = ^{l + 0{N-^))h[C], h[C]:=^<f^dx'^<f^dy^^D{x,y). (9.34) 

By taking into account all the contributions from parallel sides dx and dy (see Fig. 10), 
h[C] is calculated as, for T > i? S> 1, 

h[C] = --^^ + + A In (9-35) 

where e is the ultraviolet cutoff included to avoid the coincidence of x and y (see 
Appendix of Ref.^)). In h[C], the first term corresponds to the Coulomb potential 
in four dimensions, 

o2 1 

Vc{R) = + const. + 0{X^/N^), (9.36) 

47r K 

and the second term in h[C] corresponds to the self-energy of quark and anti-quark. 
Furthermore, if we take into account the 0{g^) correction, the coupling constant 
begins to run and the bare couplineo in (9.36) is replaced by the running coupling 
constant g = g{fi) (see, e.g., KogutEHP).0 In the topologically trivial case, therefore, 
the second term f[C] cannot give a non- vanishing string tension. 

Next, we consider the topologically non-trivial case. We begin to estimate (9.29) 
for a circular Wilson loop C with diameter R (see Fig.|Tl|). If we avoid the coinciding 
case, X = y, f[C] has a contribution only when x and y are at opposite ends of a 



* The contribution up to O(A^) in the leading order diagrams (planar diagram) in the large A'' 
expansion leads to a running coupling that differsfrom that in the usual perturbative calculation in 
the coupling constant g. 
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diameter, i.e., \x — y\ = R. Therefore, D{x,y) and G{x,y) are functions of R, due 
to translational invar iance. For any x,y£C,\x — y\ = R and 

£dx'-£*"i5(ii) = ^ = i (9.37) 



Hence we obtain 



f[C] = ^Gcix,y). (9.38) 



It is clear that f[C] is not sufficient to give a non- vanishing string tension, since 
Gc{x,y) exhibits exponential decay for large R := \x — y\. Note that Pab{x) and 
§qU = Jgdu are U(l) gauge invariant quantities, so that Gc{x,y) is also U(l) 
gauge invariant. In the large expansion, we can give a more precise estimate of 
the second term (see Appendix E). 

Finally, we consider the topologically nontrivial case of a rectangular Wilson 
loop C with side lengths R and T (see Fig. 10). In this case, we cannot give a precise 
estimate of the second term, since wc cannot perform the integration exactly. To 
leading order in the expansion, it turns out that 

Gcix,y) = G{x,y) := 2{T\biT\d{Pabix)Pcdiy))cP^N-i, (9.39) 

and that G{x, y) decays exponentially for sufficiently large \x — y\ (see Appendix E). 
Note that x and y are located on the opposite sides of the rectangular Wilson loop. 
Therefore, there exists an uniform upper bound, 

\G{x,y)\<G{R). (9.40) 

Hence, there exists an upper bound on f[C] for a sufficiently large Wilson loop such 
that T > i? > 1: 

\f[C]\ < ^G{R)h[C]. (9.41) 

Here h[C] is calculated in the same way as in (9.35) by taking into account all the 
contributions from parallel sides dx and dy (see Fig. 10). Therefore the second term 
f[C] cannot give non- vanishing string tension in the topologically nontrivial case. 

Thus, within this reformulation, the area law of the non-Abelian Wilson loop and 
the linear static potential in four-dimensional SU{N) Yang-Mills theory is realized 

if and only if the diagonal Wilson loop ( exp (i §(-, lv) \ in the two-dimensional 



NLS model obeys the area law 

( i ") )cP-- " ( Is ''^) )cPr^ ^ exp(-aoTi?). (9.42) 

In other words, the area law or the linear potential between the fundamental quark 
and anti-quark is obtained from the topological TQFT4 piece alone: 

/ exp (a <j) iojj =/ exp (i (f ui]) ~ exp(—aoTR). (9.43) 
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In any case, the derivation of the area law is reduced to a two-dimensional problem. 
It should be remarked that only the total static potential, 

V{R) = aoR - ^4^^ + const., (9.44) 
47r K 

is gauge invariant. Thus the linear potential piece alone is not gauge invariant. 
However, in the large R limit, R oo, the linear potential is dominant in V{R) so 
that the linear potential piece becomes essentially gauge invariant. 

9.2. Area law to the leading order in the large N expansion 

By the rescaling of the field (f) in the Lagrangian (9.23), another form of the 
Lagrangian of the CP^~^ model is obtained as 

CcpN-i = 9^<A"(x)a^r (x) + A(^^-{x)d,r{x) - (x)r (x))2,(9.45) 



with the constraint 



</.t(x)</)(x) := rixWix) = ^, gl := (9-46) 



It is useful to consider the Schwinger parameterization, L§) 



where 




.2 ^ -1 



U 



u'u ■.= UaUa {a = 1, ■ ■ ■ , N - 1) . (9.48) 



\UN-1 , 



Note that there is no constraint for the variable u, since the Schwinger parameteriza- 
tion automatically satisfies the constraint (9.46). We can rewrite various quantities 
in terms of u without the constraint, e.g., 

5 + 2i^dip. (9.49) 

(l + ^utu)2 g§ 

Then the U{1) gauge field V = V^dx^ in the CP^~^ model is written as 

i ok , \ , + s i ok u^du — du^u 

2^ 2 AT (1 + |^^tu)2 

We identify the complex coordinate w in the Kahler manifold with the Schwinger 
variable as 
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This leads to 



WadWa — dWaWa 



9l 



(1 - |^utu)2 N 



{UadUa — dUaUa), 



(9.52) 



and 



1 + WaWa 



(1 + 



iN 



utu) 



(1 - |^utu)2 



(9.53) 



Then we find the following expression for uj in terms of u: 



uj :- 



1 WadWa — dWaWa 

2 1 + WaWa 



n V^du — dvL^U 



2 N 



(l + ^utu)2" 



Thus the connection one-form is given by 

^ i WadWa — dWaWa 
2 1 + WaWa 



dip = UJ — dip, 



and the Abelian curvature two-form is equal to the Kahler two-form: 



dV = duj = Qk = idapdwa A dwp, 



9af3 



A5r 



,f3 - WaWf3 
A 



A 



1 + lll-u;! 



(9.54) 

(9.55) 

(9.56) 
(9.57) 



By way of the variable u, we have found that the connection one-form uj appear- 
ing in the NAST is equal to the gauge-invariant part of V. Therefore the diagonal 
Wilson loop for CP^~^ model is equal to 



exp \ i f ijj 



CP^ 



exp [I f V 



CP^ 



(9.58) 



For the CP^ model, this reduces to V (6.51) in Ref.0) for G = SU{2). 

The expectation value (exp (i y))(^pjv-i is calculaied_in Appendix D in the 
large A'" expansion in a manner based on pioneering works. ES^'ElP'Ea) The result agrees 
with the result of Campostrini and Rossi. 1^ To leading order in the 1/A^ expansion, 
the Wilson loop obeys the area law for all non-self-intersecting loops: 



exp \i <f V 



CpN- 



exp 



6tt 

'iV 



m?,|Area(C)| 



where 



m J, = jJ^ exp 



27r" 



Thus the string tension is obtained as 



67r 



^0 = ^"^^ 



(9.59) 
(9.60) 

(9.61) 
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Here is the mass of the CP^^^ field (j). The is equal to the vacuum ex- 
pectation value {(t{x)) of the Lagrange multiplier field a from the correspondence 
(T{x)<j){x)(/){x) m'^<j){x)4>{x). In the propagator of the vector field V, a massless 
pole appears. Hence, the auxiliary vector field V becomes a dynamical gauge field, 
giving rise to a linear confining potential between (p and 4>. On the other hand, 
the Lagrange multiplier field a for the constraint (9.46) becomes massive, so that it 
does not contribute to the confining potential between (p and (p. Thus we have com- 
pleted a proof of quark confinement in four-dimensional SU{N) Yang-Mills theory 
based on the Wilson criterion to leading order in the large N expansion within our 
reformulation of the Yang- Mills theory. 

§10. Remarks 

Some remarks are in order to avoid confusion. 

10.1. Calculating gauge invariant quantity in the gauge non-invariant theory 

The Wilson loop operator W^c'[^] is invariant under the gauge transformation 
denoted by U{x), which is defined by the decomposition of the variable in (7.1) 
or (7.10), since it is a gauge-invariant quantity. Therefore, the expression of the 
Wilson loop operator given by the non-Abelian Stokes theorem derived in this paper 
is also invariant under arbitrary (i.e., infinitesimal and finite) gauge transformation. 
Therefore, the variable U{x) does not contribute to the Wilson loop operator, and 
we can write 

Wc[A\ = Wc[V], (10.1) 

which is a mathematical identity. This is indeed the situation before taking the ex- 
pectation value (Wc'[-4])yM in terms of the Yang-Mills theory. However, this seems 
to contradict the claim of this paper; the area law can be derived from the contri- 
bution of the topologically non-trivial degrees of freedom expressed by the variable 
U{x), which is described by the TQFT. In fact, we have used the formula (7.19) to 
calculate the expectation value of the Wilson loop operator in § 8. Therefore, if we 
set 

g{V, U) = g{V) = Wc[Vl h{U) = 1, (10.2) 
in the expectation value {■)ym, we would obtain the inconsistent result 

{Wc[A\)ym = {Wc[V])pym{1)tqft = (H^cMW? (10-3) 

since this implies that there is no contribution from the topological part expressed 
by U and that the Wilson loop can be calculated by the perturbative part only.cl 

This apparent contradiction can be understood as follows. Elimination of the 
[/-dependence is possible only when the gauge theory itself is formulated in a gauge 
invariant way, just as the lattice gauge theory can be written in a manifestly gauge 
invariant manner without gauge fixing. In such a case, we can perform the gauge 
transformation so that a gauge invariant quantity, e.g., the Wilson loop, has no de- 
pendence on U{x) also in the expectation value. However, our strategy is totally 



* The author would like to thank Giovanni Prosper! and a referee for pointing out this issue. 
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different from tliis case. That is to say, we calculate the gauge-invariant quantity 
using the gauge fixed (i.e, gauge non-invariant) theory, where the off-diagonal compo- 
nents and the diagonal components are subject to different gauge fixing conditions. 
Therefore, the total action Sym[A,---] = Sym[-^] + Sgf+fp[-^, ■ "] obtained by 
adding the GF+FP term Sgf+fp[-^, • • •] to the Yang-Mills action S'yMi-^] does not 
have the same form as the gauge-transformed total action Sym[V, B, C, C] obtained 
by the gauge rotation U, although 5yAf[^] = S'yM[V]. In this sense, 

{Wc[A])ym + (VFc[V])yM, (10.4) 

where 

kWc{A\)YM ■■= Zyl, jvA^--- exp{i5*?l^[^, • • -WWclA]. (10.5) 

Thus it should be remarked that the formula (7.19) holds only when the field variable 
A is separated into V and U and at the same time the respective component is iden- 
tified with the topologically trivial (perturbative) and non-trivial (non-perturbative) 
components, respectively, where the V-dependent part is calculated using perturba- 
tive Yang-Mills (pYM) theory. This is an assumption of our approach. In this sense, 
Eq. (7.19) holds under this assumption and it is not a mathematical identity. In our 
approach, an arbitrary gauge transformation is not allowed under the functional of 
the expectation value, i.e., {•)ym, once the respective component is identified with 
the relevant degrees of freedom in Yang-Mills theory. (Thus, if we wish to change 
the integration variable as in ( |10.2| ), the total Yang-Mills theory should be modifed 
simultaneously. This leads to a more complicated theory.) In this paper we have 
fixed the gauge degrees of freedom for the variable U{x) using the modified MA 
gauge. In other words, we have chosen a specific field configuration of {U{x)} to 
calculate the physical quantity in such a way that the [/-dependence is controlled by 
the MA gauge fixing term. This is an essential point of our strategy. In this sense, 
the physics enters in when we consider the expectation value. 

10.2. The Gribov problem and dimensional reduction 

Dimensional reduction of the Parisi-Sourlas type was applied to theraiidQiii 
field Ising modelO and scalar .field theories with random^ternal sources o^'O' 113 
based on non-perturbative Ej^'H^) and rigorous methods. II3)'L3) In these models it 
has been recognized 0' III' that the Parisi-Sourlas correspondence (between random 
systems in d dimensions and the corresponding pure systems in d — 2 dimensions, or 
supersymmetric theory in d dimensions and the corresponding bosonic theory in d—2 
dimensions) is exact only in the case of unique solutions to the classical equation of 
motion. In fact, it was rigorously shown El^ that the three-dimensional Ising model 
in a random magnetic field exhibits long-range order at zero temperature and small 
disorder. This implies that the lower critical dimension d^ for this model is 2 (ruling 
out = 3), where the lower-critical dimension is the dimension above which [d > di) 
long-rangeferromagnetic order can exist. Therefore, this result contradicts the naive 
predictionll^P between the random field Ising model in d dimensions and the pure 
Ising model in d — 2 dimensions. 
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In gauge theories, a similar problem can in principle arise due to the existence 
of Gribov copies, although much is not known on this issue. Indeed, it is knownH' 
that Gribov copies exist for the naive MA gauge. On the other hand, it is not known 
whether the modified MA gauge with 0Sp{D\2) symmetry possesses the Gribov 
copies or not. If Gribov copies exist, first of all, the naive BRST formulation breaks 
down. Even if a modifigd BRST formulation can exist, the BRST symmetry can be 
spontaneously broken, v3 at least in the topological sector (described by the modified 
MA gauge action), where the large (or finite) gauge transformation with non-trivial 
winding number plays a crucial role. In this case, the vacuum in the topological 
sector is not annihilated by the BRST charge Qb (i-e., Qb\^)tqft / 0), whereas 
the perturbative sector (in which there exist small quantum fluctuation around an 
arbitrary but fixed topological background) is characterized by the unbroken BRST 
charge, i.e., Qb\^)pYM = 0. Even in this case, dimensional reduction should take 
place, since we have used neither the field equations nor the above property of the 
BRST charge acting on the state. In other words, the dimensional reduction is an 



origin of the spontaneous breaking of the BRST symmetry (see Refs.l7),[42D). More- 
over, the BRST symmetry can also be broken by other mechanisms, e.g., radiative 
corrections. 

In any case, the spontaneous breaking of the BRST symmetry in the topological 
sector implies that the expectation value of the gauge invariant operator depends 
on the gauge fixing parameter a for the topological sector. In this case, we can not 
conclude the a-independence of the physical (gauge-invariant) quantities. Rather, 
the value of a should be determined as a physical parameter by the theoryE^) or 
experiments. This is reasonable if we believe that the real world is described by a 
unique quantum theory and that nature is realized at an appropriate gauge fixing, 
since the classical theory does not really exist in nature, and hence ambiguities arising 
in the quantization are absent from the beginning. It is a challenge to give a definite 
answer to these questions. The author would like to thank the referee for pointing 
out the Gribov problem as an obstacle to the dimensional reduction. 



§11. Conclusion and discussion 



We have given a new version of the non-Abelian Stokes theorem for G = SU{N) 
with N > 2 which reduces to the previous resultlll' for SU{2). This version of 
the non-Abelian Wilson loop is very helpful to see the role played by the magnetic 
monopole in the calculation of the expectation value of the non-Abelian Wilson loop. 
Combining this non-Abelian Stokes theorem with the Abelian-projected effective 
gauge theory for SU{N), we have explained the Abelian dominance for the Wilson 
loop in SU{N) Yang-Mills gauge theory. For SU{N) with > 3, we must distinguish 
the maximal stability group H and the residual gauge group H, which is taken to 
be the maximal torus group H = U{1)^^^. 

In order to demonstrate the magnetic monopole dominance and the area law 
of the Wilson loop, we have used a novel refornmlation of the Yang- Mills theory 
which has been proposed by one of the authors. liSP This reformulation is based on 
the identification of the Yang-Mills theory with the perturbative deformation of a 



58 



topological quantum field theory. This framework deals with the gauge action Sym 
and the gauge-fixing action Sgf on an equal footing. The gauge action is charac- 
terized from the viewpoint of the geometry of connections. On the other hand, the 
gauge-fixing part is related to the topological invariant (Euler characteristic) deter- 
mined by a global topology. Therefore, the gauge-fixing part can have a geometric 
meaning from a global viewpoint. This point has not been emphasized in textbooks 
on quantum field theory. (See Ref.l7) for more details.) 

Our approach relies heavily on a specific gauge, the MA gauge. In spite of the 
absence of an elementary scalar field in Yang-Mills theory, the MA gauge allows 
the existence of the magnetic monopole. The magnetic monopole is considered as a 
topological soliton composed of gauge degrees of freedom, providing the composite 
scalar field. At least in this gauge, (Parisi-Sourlas) dimensional reduction occurs 
due to the supersymmetry hidden in the gauge fixing part in the MA gauge. By 
making use of the non-Abelian Stokes theorem within this reformulation of the Yang- 
Mills theory, the derivation of the area law of the non-Abelian Wilson loop in four- 
dimensional Yang-Mills theory has been reduced to the two-dimensional problem 
of calculating the expectation value of the Abelian Wilson loop in the coset G/H 
non-linear sigma model. This is the main result of this article. 

In particular, in order to demonstrate confinement of the fundamental quark in 
the four-dimensional SU{N) Yang-Mills theory in the MA gauge, we have only to 
consider the two-dimensional CP^~^ model. Prom the topological point of view, 
the Abelian Wilson loop is equivalent to the area integral (enclosed by the Wil- 
son loop) of the instanton density in the two-dimensional NLS model. This implies 
that the calculation of the magnetic monopole contribution to the Wilson loop in 
four-dimensional Yang- Mills theory was translated into that of the instanton contri- 
bution in the two-dimensional NLS model when the Wilson loop is contained in the 
two-dimensional plane, i.e. when the Wilson loop is planar. In addition, the two- 
dimensional instanton is considered as a subclass of the four-dimensional Yang-Mills 
instanton (see Ref.|l3|)). This suggests that the area law of the Wilson loop can be 
derived by taking into account the contributions of a restricted class of Yang-Mills 
instantons. A Monte Carlo simulation on a lattice will be efficient in confirming this 
dimensionally reduced picture of quark confinement. 

Moreover, naive instanton calculus (the dilute gas approximation) in the coset 
NLS model leads to the area law of the Wilson loop of the original four-dimensional 
Yang- Mills theory. Improvements of the dilute gas approximation are necessary to 
confirm the area law based on the above picture. However, this is rather difficult, as 
found more than twenty years ago.H) An advantage of the extension of the above 
strategy to SU{N) with arbitrary N is that the large N systematic expansion can 
be used in the calculation of the Wilson loop, and that the area law of the Wilson 
loop has been demonstrated. These results confirm the area law of the Wilson loop. 
It is known that the large Yang-Mills theory is related to the string theory. The 
correspondence between instanton calculus and the large expansion and a related 
issue in the large JY expansion in Yang-Mills theory will also be discussed in a 
subsequent article. e2P These considerations should shed more light on the confining 
string picture. 
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Appendix A 

Normalization of the Coherent State 



In this appendix we derive the normahzation facter N of coherent state l^^,^) 
which parametrizes G/H, 



N := {A\ exp [f„^„] exp [r^^./^] |^) 



(A-1) 



1 1 



E ITZi^"^ • • • • • • • • • ^«^-^-/3i • • • E_f^L \A) .(A-2) 



K,L=0 



A.l. SU{3) coherent state 

We use the positive root a*-*^ and Dynkin index [m, n] in SU{3), which are defined 
by Fig. 3 and ( 3.13| ). This definition may be rewritten as 



2 ' 2 



(A-3) 



For [m,n] = [m,0] (resp., [0, n]), a and /3 run over 1,2 (resp., 2,3). Then the 
corresponding coherent state parametrizes CP^. From the orthogonality of the states 
which span the representation space, the terms contributing to in (A-1) must 
satisfy the condition 



+ . . . + a-^ - 



0. 



(A-4) 



In the [m, 0] case, this condition implies that the number of a^^^ and a^^^ in K positive 
roots is equal to that of a^^^ and a*^^^ in L negative roots. Since [Ei,E2] = 0, we 
have only to estimate the terms 



Nk^i '■= {A\ El ... El E2 ■ ■ ■ E2 E—2 . . . E-2 E-i . . . E—i \ A) . 

k 



(A-5) 



I 



We begin with the term 



0,1 



{A\ E2...E2E.2--- E^2 \A) = No^i (a - la^'^^\A - I 



la 



(2)\ 



(A-6) 



where we have used in the last equality the fact that E-2 ■ ■ ■ E-2 |^) has a weight 
A — la^"^^ and is proportional to the state A — la^"^^), which is normalized as 



(^A-la^MA-la^^^^ = l. 



Exchanging the rightmost E2 with E-2 and using 

[E2,E.2]=a'^^^ -H, 



a(2) .h\A- ja^^)) = - j) \a - ja^^)) , (0 < i < m) 



(A-7) 



(A-8) 
(A-9) 
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we obtain the recursion relation for Nq^i, 

No,i = 1^ - (/ - 1)| No,i-i + {A\ E2... E2E_2E2 E_2 • • • ^-2 |^) 

/-I i-i 

+ f -G-2) + --- + f)iVo,.-i 
-iVo,/-i, 



m 

l{m-l + l) 



(A- 10) 



where we have used E2 \A) = 0. Prom this relation and from Nq^q = {A\A) = 1, we 
obtain 



0,1 



Similarly, 



N, 



k,0 



2' ' 



2k 



(A-ll) 



(A-12) 



For the general terms (A-5), we have 



Nk,i = iVfc,o (a - ka^'^ 



Eo . . . Eo E. 



2 -C/-2 • - 



E. 



A - fca(^) 



k{m-k-l + l) 



A - ka^^^ 



Eo . . . Eo E—o . . . E. 



l-i 



l-i 



2'^ 2' 

Finally, the normalization factor N is given by 



A-ka'^^y 

(A-13) 



1 



k + r 



T2 . . . T2T1 . ^ . ri 7-2 . T2T1 .^.Tl Nk^l 
I k I k 



n n—k 

EE( 

k=0 1=0 



T2T2 \k ( ^\ (m- k 



k 



I 



= (1 + (wiwi) + {w2W2)Y^ = exp Kcp2{w,w), 



(A-14) 



where wi := W2 := T2/V2, and K(jp2{w) is the Kahler potential of CP^. 

For F2, i-e., mn / 0, a and /3 in (A-1) run over 1, 2, 3. Here we apply almost the 
same steps as in the case of CP^. The result is 

3 3 

N := {A\ exp ^ [fe^E^] exp ^ [r^E.p] \A) 

a=l 13=1 
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\/2 V2 4 ' y V V2 V2 4 

= (l + l^il +\W2\ )) (l + |lf3| +\W2-WiW3\ 1 

where := ti/\/2, W2 := T2/\/2 + rir3/4, tt;3 := T-^j^f^. 

A. 2. SU{N) coherent state 

In SU{N), the simple root a* and the Dynkin index [mi, • • • rriN-i] are defined 
by (3.47) and a* • yl := mi/2. In general, the roota_E'±Q belonging to the coset group 
G/H are not orthogonal to the highest weight Axi We restrict our consideration to 
only the case of [m,0, ••• ,0], in which case only is not orthogonal to A. Any 
positive root a is given by a linear combination of simple roots a*. The positive 
roots Oi in (A-1) are not orthogonal to A, so they must include a^. It turns out that 
such roots are given by cij := + • • • + a* (i = 1, • • • — 1) , which are realized by 
{Eajj^i = 6i^kSi+i,i/ in the N representation. The corresponding coherent state 
parametrizes CP^~^. We can easily extend the CP^ case to that of CP^~^ using 
the fact that [£"0,^ , £'5^. ] = and on ■ A = m/2. (This situation is same as in the 
5C/(3) case, i.e., [Ei,E2] = 0, and a^^) • /= a^^) . A = m/2.) Thus we obtain 

N-l N-l 

N := {A\ exp ^ [f.^^J exp ^ [r.E.o.,] \A) 

i=i 

/ N-l 
1^1 + ^ {WiWi)^ 

exp Kqpn-i{'w,w). (A- 16) 

Appendix B 

From CP^ to CP^ 




For an arbitrary element ((/>!, (j)2, (ps) of W, 

= - (0,0,0), (B-1) 

the entire set of ratios (pi : (f)2 ■ 03 is called the complex projective plane and is 
denoted by P'^{C) or CP"^: 

{a(pi : a4>2 ■ "(As) = {<Pi '■ (p2 ■ (ps), a e C, a ^0. (B-2) 

Defining the subset Ua{a = 1, 2, 3) of CP^ as 

Ua = {(01 : 02 : 03) e CP2; 0, / 0} C CP^ (B-3) 

we observe that 

(01 : 02 : 03) ■■ 



1 : ^ : ^ ) G f/i. (B.4) 



* From the definition of the coherent state (2.9), we have E±a \A) 7^ and [Ea,E-a] \A) 
a-A\A)j^O. 
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The mapping ipi from Ui to C^, 

'PI : {<t>i : h : 4>3) eUi^(^:^)e (B-5) 

V<Pi 91/ 

is a bijection, i.e., a surjection (onto-mapping) and injection (one-to-one mapping). 
The inverse mapping is given by 

ip^^:ix,y)£C^^{l:x:y)GUi. (B-6) 

Similarly, the following maps (p2 and (ps are also bijections from U2 and Us to C^: 

f2 : (01 : 02 : eU2^(^:p-)e (B-7) 



V'S : (01 : 02 : 03) e C/3 ^ ^ : ^ G C^. (B-8) 

V 03 03 / 

Since 

CP2 = [/, u {(0 : 02 : 03)}, (02 : 03) 7^ 0, (B-9) 

(02 : 03) determines a point in CP^. Conversely, for a point (61 : 62) in CP"*", 
(0 : 61 : 62) defines a point in CP^ — Ui. Thus the map 

(0 : 02 : 03) e CP' -Ui^ (02 : 03) G CP^ (B-10) 

is one-to-one and onto. Then we can identify CP'-C/i with CP^ as CP'-J/i = CP^ 
On the other hand, the identification of Ui and C^, Ui = C^, by the mapping 
cpi leads to 

CP' = ?7i U CP^ ^ C' U CP^ (B-11) 

Using CPi = C^U{{0 : 1)}, we can write CP' = C'uCiu{(0 : 1)} = C'UC^UCO. 
CP' — C/i is called the "line at infinity" and is denoted by £00, CP' — Ui = loo- 
Thus we can also write 

CP' = ?7iU£oo, C/i = C', ^oo = CP\ (B-12) 

where C' is called the "complex affine plane". The homogeneous coordinates (0i : 
02 : 03) are related to an element (x,y) in the affine plane as 

^ = ^' y = ^^ ^00 = {01 = 0}. (B.13) 

01 01 

CP' is obtained as CP' = W U V U W by gluing the three affine planes U, V and 
W using the relations 

^ V2 wi I I Ul 
Ul = — , U2 = — , ^^1 = , V2 = , Wi = , W2 = , (B-14) 

Vi Vi W2 W2 U2 U2 

where {ui,U2), {vi,V2) and (^1,^2) are the coordinates ofU,V and W defined in 
terms of homogeneous coordinates of CP' as 

^"^'"^^=(1'!)' ^^^'^^^=(1'!)' ^"^'"^^=(l^'l)- ^^-''^ 
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Note that the method of gluing affine planes is not unique. In fact, CP^ can be 
obtained from the following gluing: 

1 1 Uo Vl 1 



W2 i L U2 Vi 

Ui = , U2= , Vi = , V2 = , Wi = —, W2 = —. 

Wi Wi Ui Ui V2 V2 

Appendix C 

Nonlinear Sigma Model of the Flag Space 



C.l. Correspondence between SU{N)/T and SL{N,C)/B 

When an element ^ in F/v_i is expressed by the complex coordinate, 



(B-16) 



/ 1 Wi W2 

1 Wn+l ■■■ 
1 W2n 





Vo 



Wn \ 
W2n-l 
W3n-3 



1 Wn{n+l)/2 

1 



(C-l) 



we find det^ = 1. Hence ^ is an element of SL{N,C). There is the isomorphism 

SU{N)/T = SL{N,C)/B. However, ^ in this form is not necessarily unitary. The 
corresponding unitary matrix V G SU{N) is obtained as follows. First, ^ as an 
element of SL(N, C) is expressed in terms of the column vectors: 

i = {Ei,E2,---,EM)^ SL{N, C) = SU{Nf. (C-2) 

By applying the Gramm-Schmidt orthogonalization, we can obtain a set of mutually 
orthogonal vectors {E[, £"2, • • • , E'j^) from (Ei,E2, • • • , Ej\f) as 



E'l :— El, 
E'2 ■= E2 - 



{E2,E[ 



iE[,E[ 



p/ ._ p iEN,Ej^_^) {En, El) , 

l-fc'iV-i'-^Ar-iJ \Ei,E-^) 



(C-3) 



where the inner product is defined by {Ei,Ej) := EJ ■ Ej. Using the normalized 
vectors ej := E'-/\\E'-\\, we obtain an element in SU{N), 



V={ei,e2,---,eN)eSU{N). 
In fact, det V = 1 and V is unitary, since 



(C-4) 



(ei,e2,- ■ ■ ,eAr), (V^V)i 



-T 

e^ ' c-j 



(Ci , Cj 



(C-5) 



'AT. 
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In general, the unitary matrix V is related to its complexfication V'-^ by 

V = V^B, V G SU{N), e SL{N, C), (C-6) 

where B is an upper triangular matrix. This is nothing but the Iwasawa decomposition.^ 
Since the upper triangular matrices form a group, we have 



(C-8) 



where B' = B~^ is also upper triangular. This implies ^ = VB by the above 
construction. Therefore, V is indeed the element of SU{N) corresponding to ^. 
Note that the multiplication by the matrix B leaves the highest- weight state 
invariant, so that 



^\A) = VB\A) = V\A). 
For example, d = E[/\\E[\\ = Si/||^i||, i.e.. 



(C-9) 



1 



/ 1 \ 

— Wl 



ei 



Z\V2 



N 



A:=l + WwW"^ = 1 + \^^\ 



(C-10) 



\-wnJ 
The Mauer-Cartan form is 



a=l 



-T 

ei 



f Cidei eide2 
°2dei e2de2 



-T 



\ejrdei eLde2 



ei dcN \ 
eJdcN 

eJrdeN J 



which can be decomposed into diagonal and off-diagonal parts as 



N 



V ^dV = Y, efdeJii + ^ e^debE, 



ab- 



1=1 



In the fundamental representation, we have 



= {A\iV-UV\A) = i{V-'dV)ii = ie{ de 

= {A\V-^T^V\A) = ef (r^)ei. 
Hence, substituting (C-10) into (CT3), we obtain 

i WadWa - dWaWa 



A 



* Any element Qc £ G may be factorized as 

gc = gh, geG, beB 
in a unique fashion, up to torus elements that are common to G and B. 



(C-11) 
(C-12) 

(C-13) 
(C-14) 

(C-15) 
(C-7) 
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The Lagrangian density of the coset G/H, i.e. the flag NLS model, 

Cnlsm = ^ivG/H{iV-^d^,ViV-^d^V), (C-16) 
can be written in terms of the off-diagonal elements as 

^NLSM = Pg Yl i^t^)ab{^iJ.)ab = Pg Yl i^a, dfj,eb){ea, dfj,eb), (C-17) 
a,b:a<b a,b:a<b 

where wc have used (rj, e.j) := ej ■ ej = 5ij. 

Especially, the Lagrangian of the CP^~^ model is obtained as a special case of 
(C-17) as follows. Using the definition (4.37), 

= {UT^U^)n = Uia{T%bUib = U^a{T^)abUl, (C-18) 
and C/i/t = 1, we find 

N N 

= Y{iUd^U^)uiiUd^U^)u = ^(ei, 5^66)2. (C-19) 

6=2 b=2 

If we use (4.36) ((^a = f^ia = Ul^), we obtain another expression, 

5^n^5^n^ = di,4> ■ d^^ + (0 • d^^){^ ■ 5^0), (C-20) 
where we have used • ^ = 1. Thus the Lagrangian of CP^~^ model is obtained as 

CcpN-, = ^d^,n-d^T, (/x=l,---,d) (C-21) 



2 



N 



'-Yieud^ebf (C-22) 



2 



6=2 



where 



9a0) ■= - <Pc>4>l3- (C-24) 

This agrees with the Lagrangian obtained from the Kahler potential, 

-Ccp--! = ^ga0Hdf.w"d^w^ (C-25) 

9ap{w) = , Mi„..m9A9 -■ (C-26) 



2 -a/3 

with 



(l + |||u;|||2)2 

The explicit construction of V is given in the following section 
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C.2. SU(2) 

For G = SU{2), wc have 

It is easy to see that 

{Ei,Ei) = l + ww, {E2,Ei) = -w. (C-28) 

Hence we obtain 

ei = ^-^/^(_'^), e2 = zl-^/^(f), Z^:=l + H^ (C-29) 



and 



The elements of the one-form V ^dV are 

ejdei = -A''^(wdw — wdw), (C-31) 

eldei = -A-'^dw. (C-32) 
The Lagrangian of the Fi = CP^ model reads 

jO-nlsm = -(igA~'^dfj_wdf,w = -(3g- ^ d^wd^w. (C-33) 

(1 + |w| ) 

C.3. SU(3) 

For G = 5'C/(3), we have 

1 0^ 

-wi 1 I , (C-34) 



-W2 W3 1 , 



or, 



Ei=\ -wi\, ^2 = ( 1 J , ^3 = ( j . (C-35) 

-W2/ \W3/ \ 1 / 

Then we have 

{Ei,Ei) = I + W1W1 + W2W2, {E2,Ei) = —wi — wsW2,- ■ ■ ■ (C-36) 
A straightforward calculation leads to F = (61,62,63) G SU{3), with 

ei = (zli)-V2 

\-'W2. 
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(iVl + WsW2 
1 + |ii;2p - W1W2W3 
— W2W1 + W3 + ws\wi\'^ ^ 

e3 = {A2)-^'^\ -ws j, (C-37) 

where 

Ai := 1 + \wif + \W2\'^ , A2 := I + \W2 - wiwsf + Imsl"^ . (C-38) 
The off-diagonal elements of the one-form V~^dV are 

e^dei = (zAi)"^(Z\2)"^/^[(l + \w2\'^ - wiW2W3)diDi + {-W2W1 +W3 + W3\wi\'^)dw2]., 
ejdei = iAi)-^/\A2)-^/^[dw2 - iU3dwi], 

elde2 = {Ai)-^/^{A2)-^[{wi + wsW2){dw2 - w^dw-C) - A-^^dw^]. (C-39) 



For CP^, we have 



ei = (Z\i)-i/2 j I , (C.40) 

e2 = {A^A2)-^'^ ( 1 + J , (C-41) 

(W2\ 

63 = (/A2)-^/' ( 1 , (C-42) 

where 

Z\i := l + |w;i|2 + |u;2|2 = l + A2 := I + \w2\^ . (C-43) 

The off-diagonal elements of the one-form V~^dV are 

e^dei = (zli)~^(Z\2)~^/^[(l + \w2\^)dwi - W2Widiij2], 

e^dei = {Ai)-^/HA2)-'/^[dw2], 

ejde2 = {Ai)-^/\A2)-^[widw2]. (C-44) 
The Lagrangian of CP^ model is given by 

'Ccp^ = f EK^i'^Meo)!' (C-45) 
^ 6=2 

= Y (i^||[^|||2)2 [(l + \w2\^)^^.Wl^^Wl + (1 + kl|')5^«;25^lt)2 

-iiJiW2dij,widf^ijU2 - wiW2dij,W2dfj,iBi]. (C-46) 



Appendix D 

Large N Expansion of the CP^"^ Model ■ 



The generating function of the CP^ ^ model is defined by 



Z[J,J,J^] := J Vml[5 (^\(Pix)\' 



N 

7o 



X exp j-S" + J (fx[J ■(f) + (j)-J + Jj^F/,]| , 
where S is the action of the CP^~^ model, 



(D-1) 



(fx 



d^x 



% 



(D-2) 
(D-3) 



and the auxiliary vector field is defined by 



Introducing the Lagrange multiplier fields a{x) and An{x), we can rewrite as 



IlMl'^C 



7o 



= Va P^^iCxp 



exp I j d^x 



^V^{x)V^{x) + J^{x)V^{x) 



}• 



9l 



+^A^[i{ct^.d,<l^) + J,]-f^J^J, 



(D.5) 



where we have inserted the mass term m?\(l)\'^ for later convenience and have chosen 
a specific normalization for the field a. Then we obtain 



Z[J,J,J^] := j V(f)V(j} Jt^(^ J VAfj^expl - J 



+ / d^x 



X 



9l 



■ ^B(p-\ 

9o 



N 



(D-6) 



where 



Ab := -D^D^ + rr?- := + -^A^{ 



N 



(D-7) 
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This theory has global SU (N) invariance corresponding to rotations of (pa- Moreover, 
it has local U{1) gauge invariance under the transformation 



<PUx) = e'^^^Ua{x), {a = l,.--,N) 
A'^ix) = A^{x) - VNd^,A{x), 
u'{x) = (t[x). 

We can perform the integration over (j) and to obtain 



+ 



X 



1 (? ' 

7 1 A T _ 7 7 



where 



Seff := NTrln Ab + / d^xa{x). 

The effective action can be expanded in a power series of X/N: 

oo 



(D-D) 



(D-10) 



(D-ll) 



n=l 



The diagramatic representation of this expansion is given in Fig. 13 using the rule 
given in Fig. 12. 

First, the order iV^/^ term corresponds to the diagrams (a) and (b) in Fig. 13: 



9o J 

(fq 



= ia{0) 



1 r d^q 1 
^"i (27r)2 g2 + ^2 



where we have used the Fourier transformation 

a{p) = J d'^xe-'P'^aix). 



(D.12) 
(D-13) 

(D-14) 



The integral in (D-13) is ultraviolet divergent. It can be regularized by introducing 
the cutoff A. The saddle point condition S^^^ = requires the bare coupling constant 
go to vary with the cutoff A according to 



1 



(D-15) 



g^{A) Att VP?' 

In other words, if the bare coupling varies with respect to the cutoff according to 

1 



1 \ 

= — In ■ 



(D.16) 
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the divergences cancel in S^^\ This imphes the asymptotic freedom of the CP^ ^ 
model. Imposing the condition S^^^ = 0, therefore, we obtain 



= exp 



47r 



(D-17) 



Next, the order term corresponds to the diagrams (c),(d) and (e) in Fig. 13: 
^^'^ = ^/'^'^Z d^y[<y{x)r{x,y)a{y) + A^{x)r^,{x,y)A,{y)], (D-18) 

where the Fourier transformation of r{x,y) and r^j^i,{x,y) are respectively given by 
(see Fig. 13) 



r{p) 



:ln 



_|_ 4j^2 _|_ 



and 



(D-19) 



(27r)2 g2 + 



p2 



ip' + 4m")r(p) - - 



(27r)2 (p2 + m2)((p + g)2 + m? 
1 



(D-20) 



In the neighbourhood of = 0, we have 



r{p) = 



1 



247rm4 



127rm2 



V- p2 

(D-21) 



Thus we obtain the low-energy effective action, 

(D-22) 

It is important to note that the kinetic terms for the Lagrangian multiplier fields are 
generated. The field a becomes massive, while the field ^//(= VNVfj,) is massless. 
The line integral in the Wilson loop can be rewritten as 

£ d^MO = J d^zV^{z)J^{z) =: {V^, J^), (D-23) 

1 (zeSc) 



where 



Mz) = i d^^,5\z -0 = €^,d,^{z), #(z) 



c 



{ziSc) 



(D-24) 
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Here Sc is the area bounded by the loop C. Up to leading order, we can perform 
the Gaussian integration to obtain 



exp 



C 



CP^ 



[0,0,0] J VaVA^e-^^^^ exp 
2 



const exp 



N 



Thus we obtain the area law. 



exp 



CpN-1 

since J = and A := d5 + 6d, and hence 



const exp 



67rm 



\Sc\ 



(D-25) 



(D-26) 



(J^,Z\-ij^) := J d\ J d^yJ^{x)A-\x,y)J^{y) = = \Sc\. (D-27) 

The dynamically generated gauge field produces a long-range force with a linear 
potential that confines the (ps. Both global SU{N) and local U{1) symmetries are 
unbroken in two dimensions due to the Coleman theorem. In dimensions D > 2, it 
has been shown LJ) that there is a critical point gc such that for g < Qc SU{N) and 
[/(I) symmetries are broken, 7^ 0, while {a) = 0. In this phase (/> is regarded 
as the Nambu-Goldstone particle, since = 0. A massless vector pole exists in 
the propagator {A^{x)Ay{{))) .V3 For g > gc-, on ther other hand, the SU{N) and 
U{1) symmetries are exact, implying {(j)a) = 0, and (a) / 0. In this phase, cj) is 
massive, m,^ = {a). For D = 2, gc = 0. For more details on large results, see 
Refs]6|),|^,|8|),|6|). 

Appendix E 

Large N Estimation 



If we write n^{x) in terms of the CP^ ^ variable 
we obtain 

= Pah{x)Pcd{y){T^)ab{T%d 



UE as 



(E-l) 



(0(x) - ^{y)){m - ^{y)) - -(0(x) - ^{x)){m - 0(y)) 



where we have used 



E 

A=l 



^adhc 



N 



^abScd ■ 



(E-2) 



(E-3) 
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Hence we obtain 

1 



n{x) ■ Ta.{x) = n'^{x)n"^{x) = ^ 



(4>(x) ■ 4.(x)r. (E-4) 



The contraint <p{x) ■ (^{x) = 1 leads to n(a:;) • n(a;) = ^ ^ ~ 7! ■ ^he expectation 
value reads 

2{n\x)n\y)) = • </.(y))(0(x) • ^{y))) - l{\cl>{x)\^\^{y)f). (£-5) 

The factorization in the large A'^ expansion leads to 

2{n^{x)n^{y)) - • </>(2/))(</>(x) • ^(y)) - ^(|</>(x)p)(|0(y)|2), 

= \{4>{x) ■ m)]' - l(|</,(x)p)(|0(y)p). (E-6) 

The large expansion shows that the field (p becomes massive, so that the two-point 
function {4>{x) ■ 4>{y)) exhibits exponential decay. 

To leading order in the large N expansion, the correlation function Gc{x,y) 
defined by (9.32) in the CP^^-'^ model reads 

Gc{x,y) = G{x,y) := (^Sadhc - j^Kb^cd^ {Pab{x)Pba{y)) cp^-^- (E-7) 
Note that P is a projection operator with the properties 

p2(x) = P{x), trP(a;) = 1. (E-8) 
The composite operator P is gauge invariant. SU{N) invariance implies that 

{P^{x)) = {4>a{x)Ux)) = ^6ab. (E-9) 

It is easy to show that 

G{x,y) = {(^a{x)(l)bix)4>biy)<l>aiy)) - 

= {D^{x,y)) + j^{D{x,x)D{y,y)) - ^ 



where {D{x,x)) = 1. 



{D\x,y)) + ^{D{x,x)D{y,y))conn, (E-10) 



Acknowledgments 



The aTithors woTild like to thank MinorTi Hirayama for sending valuable com- 
ments on the non-Abclian Stokes theorem derived in this paper after submitting this 
article for publication. This work is supported in part by a Grant-in-Aid for Scien- 
tific Research from the Ministry of Education, Science and Culture (No. 10640249). 



73 



References 



[1] 
[2] 

[3] 



[4] 
[5] 

[6] 



[7] 

[8] 
[9] 



[10] 



[11] 
[12] 
[13] 
[14] 
[15] 
[16] 
[17] 
[18] 

[19] 
[20] 



M. Gell-Mann, Phys. Lett. 8 (1964), 214. 
G. Zweig, CERN Rep. 8419/TH 412. 

C. N. Yang and R.L. Mills, Phys. Rev. 96 (1954), 191. 
R. Utiyama, Phys. Rev. 101 (1956), 1597. 

R. Shaw, The problem of particle types and other contributions to the theory of elementary 
particles, Cambridge University Ph.D. thesis (unpublished). 

G. 't Hooft, in Conference on lagrangian field theory, Marseille, 1972 (unpublished). 

D. J. Gross and F. Wilczek, Phys. Rev. Lett. 30 (1973), 1343. 

H. D. Politzer, Phys. Rev. Lett. 30 (1973), 1346. 
K. Wilson, Phys. Rev. DIO (1974), 2445. 

Creutz, Phys. Rev. Lett. 45 (1980), 313 
Creutz and K.J.M. Moriarty, Phys. Rev. D26 (1982), 2166. 
Nambu, Phys. Rev. DIG (1974), 4262. 

't Hooft, in High Energy Physics, ed. by A. Zichichi (Editorice Compositori, Bologna, 
1975). 

S. Mandelstam, Phys. Report 23 (1976), 245. 

A.M. Polyakov, Phys. Lett. B59 (1975), 82; Nucl. Phys. B120 (1977), 429; Gauge fields 
and strings (Harwood Aca demic Publishers . London, 1987). 
N. Seiberg and E. Witten. |iop-th/9407"087| , Nucl. Phys. B426 (1994), 19. 
T. Suzuki, |hep-lat /9506016]" 



R.W. Havrnakc r, |hcp-lat/9510035 
A. DiGiacomo, |hcp-th/9603029|. 



hep-th/9710205 



M.I. Polikarpov, ^iep-lat/960"902C . 
M.N. Che rnodub and M.I. Pohkarpov. 
G.S. Bah, ^iep-ph/980935l[ 
G. 't Hooft, Nucl. Phys. B190 (1981), 455. 

A. Kronfeld, G. Schierholz and U.-J. Wiese, Nucl. Phys. B293 (1987), 461. 
A. Kronfeld, M. Laursen, G. Schierholz and U.-J. Wiese, Phys. Lett. B198 (1987), 516. 
K.-I. Kondo and Y. Taira, |hep-th/9906129| . Mod. Phys. Lett. 



A13 (2000), 366. 



K.-I. Kondo, icp-th/970910c , Phys. Rev. D57 (1998), 7467 
K.-I. Kondo, iep-th/980102^ , Phys. Rev. D58 (1998), 105019. 
K.-I. Kondo, iep-th/980313o , Phys. Rev. D58 (1998), 085013. 
K.-I. Kondo, lep-th 798051^ , Phys. Rev. D58 (1998), 105016. 
K.-I. Kondo, iep-th/981016'/ , Phys. Lett. B455 (1999), 251. 
K.-I. Kondo, iep-th/990404£ , Intern. J. Mod. Phys. A , to be published 
M. Faber, J. Greensite and S. Olejnik, ^iep-lat/980700^ ; |hep-lat/9710039| , Phys. Rev. D57 
(1998), 2603. 

R. Wensley, Poster at Lattice 99, Pisa, Italy. 
M.B. Halpern, Phys. Rev. D19 (1979), 517. 
N.E. Bralic, Phys. Rev. D22 (1980), 3090. 
I.Ya. Aref'eva, Theor. Math. Phys. 43 (1980), 353 

Yu.A. Simonov, Yad. Fiz. 50 (1989), 213 [Sov. J. Nucl. Phys. 50 (1989), 134]. 
B. Broda, Phys. Lett. B262 (1991), 288. 
B. Broda, J. Math. Phvs. 33 (1992), 1511. 



F.A. Lunev, 
M. Hirayama and S. 



hep-th/9609166t Nu cl. Phys. B494 ( 1997) 



Matsubara, 



liep-th/9712120, Prog, 



433. 
Theor. 



Phvs. 



99 (1998), 691; 
hep-th/9806098|. Prog. Theor. Phys. 



J.S. Rno, [hep-th/9903"22lt ^iep-th/9910"l7^ , J. Math. Phys. 40 



M. Hirayama, M. Kanno, M. Ueno and H. Yamakosh 
100 (1998), 817. 
R.L. Karp, F. Mansouri 
(1999), 6033. 

[21] D.I. Diakonov and V.Yu. Petrov, Phys. Lett. B224 (1989), 131; tiep-th/9606104 
[22] D.H. Feng, R. Gilmore and W.-M. Zhang, Rev. Mod. Phys. 62 (1990), 867. 
[23] A.M. Perelomov, Phys. Report 146 (1987), 135. 

A.M. Perelomov and M C. Prati. Nucl. Phvs. B258 (1985), 647. 
[24] M.-H. Kim and P. Oh, ^iep-th/9504001 . 
[25] Z.F. Ezawa and A. Iwasaki, Phys. Rev. D25 (1982), 2681. 
[26] T. Suzuki and I. Yotsuyanagi, Phys. Rev. D42 (1990), 4257. 



74 



[27] G. 't Hooft, Nucl. Phys. B79 (1974), 276. 

A.M. Polyakov, JETP Lett. 20 (1974), 194. 
[28] L.G. Yaffe, Rev. Mod. Phys. 54 (1982), 407. 

[29] S. Coleman, 1/N in Aspect of Symmetry (Cambridge Univ. Press, New York, 1985) . 
[30] S.R. Das, Rev. Mod. Phys. 59 (1987), 235. 

[31] P. Rossi, M. Campostrini and E. Vicari, Phys. Reports 302 (1998), 143. 
[32] G. 't Hooft, Nucl. Phys. B72 (1974), 461. 

[33] K. Ueno, Introduction to algebraic geometry (in Japanese), (Iwanami, Tokyo, 1995). 
[34] M.A.B. Beg and H. Ruegg, J. Math. Phys. 6 (1965), 677. 

E. Chacon and M. Moshinsky, Phys. Lett. 23 (1966), 567. 

T.J. Nelson, J. Math. Phys. 8 (1967), 857. 

D.F. Holland, J. Math. Phys. 10 (1969), 531. 
[35] M. Nakahara, Geometry, Topology and Physics, (lOP, Bristol, 1990). 

[36] M.B. Green, J.H. Schwarz and E. Witten, Sueprstring theory. Vol. 2: loop amplitudes, 
anomalies & phenomenology (Cambridge Univ. Press, Cambridge, 1987). 

[37] H. Georgi, Lie algebras in particle physics, from isospin to unified theories (Ben- 
jamin/Cummings, London, 1982). 

[38] R.F. Picken, J. Math. Phys. 31 (1990), 616. 

[39] D.P. Zelobenko, Compact Lie groups and their representations (AMS, Providence, 1973). 
[40] A. Perelomov, Generalized coherent states and their applications (Springer, New York, 
1986). 

[41] F. Wilczek and A. Zee, Phys. Rev. Lett. 52 fl9 84\ 2111. 

[42] K.-I. Kondo and T. Shinohara, [hep-th/000415^ , Phys. Lett. B491 (2000), 263 
[43] K.-I. Kondo, in preparation 



[44] L. Faddeey and A.J. Niemi, [hep-th/9807069| ,; Phys. Rev. Lett. 82 (1999), 1624; [hep- 

th/981209C| , Phys. Lett. B449 (1999), 214. 
[45] P.A.M. Dirac, Proc. Roy. Soc. London A133 (1931), 60; Phys. Rev. 74 (1948), 817. 
[46] T.T. Wu and C.N. Yang, P hys. Rev. D12 fl 975'). 3845; Nucl. Phys. B107 (1976), 365. 
[47] M. Hirayama and M. Ueno, |hep-th/9907063t Prog. Theor. Phys. 103 (2000), 151. 
[48] E.J. Weinberg, Massless and massive monopoles carrying nonabelian magnetic charges. 



hep-th/9610065 . 
[49] A. Sinha, Phys. Rev. D14 (20 16). 1976 

[50] M. Quandt and H. Reinhardt, [hep-th/9707"l85| . Int. J. Mod. Phys. A13 (1998), 4049. 
[51] K.-I. Kondo and T. Shinohara, Renormalizable Abelian-pr oiected effective gauge theory 

derived from Quantum Chromodynamics, CHIBA-EP-121, hep-th/0005125| , and paper in 

preparation. 

[52] K. Amemiya and H. Suganuma, Phys. Rev. D60 (1999), 114509. 
[53] M. Liischer, Phys. Lett. B78 (1978), 465. 
[54] D.J. Gross, Nucl. Phys. B132 (1978), 439. 

[55] B. Berg and M. Liischer, Commun. Math. Phys. 69 (1979), 57. 

V.A. Fateev, I.V. Frolov and A.S. Schwarz, Nucl. Phys. B154 (1979), 1; Sov. J. Nucl. 
Phys. 30 (1979), 590. 

A. P. Bukhvostov and L.N. Lipatov, Nucl. Phys. B180 (1981), 116 

P.G. Silvestrov, Sov. J. Nucl. Phys. 51 (1990), 1121. 
[56] A. D'adda, P. di Vecchia and M. Liischer, Nucl. Phys. B146 (1978), 63. 
[57] A. D'adda, P. di Vecchia and M. Liischer, Nucl. Phys. B152 (1979), 125. 
[58] E. Witten, Nucl. Phys. B149 (1979), 285. 

[59] G. Miinster, Phys. Lett. B118 (1982), 380. Nucl. Phys. B218 (1983), 1. 
[60] M. Campostrini and P. Rossi, Phys. Rev. D45 (1992), 618; Phys. Lett. B272 (1991), 305. 
[61] D. Gross, ^iep-th/9212"l4E| , Nucl Phvs. B400 (1 993), 161 



D. Gross a nd W. Taylor IV, ^ep-th/9301068t Nucl. Phys. B400 (1993), 181; [hep- 
th/930304e| , Nucl. Phys. B403 (1993), 395. 
[62] J.B. Kogul!, Rev. Mod. Phys. 55 (1983), 775. 
[63] J. Schwinger, Phys. Lett. B24 (1967), 473, Phys. Rev. 167 (1968), 1432 
[64] K.-I. Kondo an d Y. Taira. Nucl . Phys. Proc. Suppl. 83, 497 (2000). 

[65] A.M. Polyakov, |hep-th/9607049|. Nucl. Phys. B486 (1997), 23; hep-th/971002, Nucl. Phys. 

B Proc. Suppl. 68 (1998), 1; [hep-th/9809057| . Int. J. Mod. Phys. A14 (1999), 645. 
[66] M. Stone, Nucl. Phys. B152 (1979), 97. 



75 



[67] H.E. Haber, I. HinchlifTe and E. Rabinovici, Nucl. Phys. B172 (1980), 458. 
[68] A.J. Macfarlane, Nucl. Phys. B152 (1979), 145. 

[69] W.A. Bardeen, B.W. Lee and R.E. Shrock, Phys. Rev. D14 (1976), 985. 
[70] I.Ya. Arefeva and S.I. Azakov, Nucl. Phys. B162 (1980), 298. 
[71] T. Kugo and P. Townsend, Phys. Lett. BlOl (1981), 82. 

S. Weinberg and E. Witten, Phys. Lett. B96 (1980), 59. 

T. Kugo, Phys. Lett. B109 (1982), 205. 

T. Kugo, H. Terao and S. Uehara, Prog. Theor. Phys. Suppl. 85 (1985), 122. 

T. Kugo, Soryushiron Kenkyu (Kyoto) 71, E78-90 (1985); Prog. Theor. Phys. 66 (1981), 

2249. 

M. Bando, T. Kugo and K. Yamawaki, Phys. Reports 164 (1988), 217. 
[72] Y. Imry and S.K. Ma, Phys. Rev. Lett. 35 (1975), 1399. 
[73] J.Z. Imbrie, Commun. Math. Phys. 98 (1985), 145. 

J. Bricmont and A. Kupiainen, Phys. Rev. Lett. 59 (1987), 1829. 
[74] J. Cardy, Phys. Lett. B125 (1983), 470. 
[75] A. Klein and J. Perez, Phys. Lett. B125 (1983), 473. 
[76] A. Klein, L. Landau and J. Perez, Commun. Math. Phys. 94 (1984), 459. 
[77] G. Parisi and N. Sourlas, Nucl. Phys. B206 (1982), 321 

[78] J. Zinn- Justin, Chapter 16: Slavnov- Taylor and BRS symmetry. Stochastic field equations, 

in Quantum Field Theory and Critical Phenomena (Oxford Univ. Press, 1993). 
[79] K. Fujikawa, Nucl. Phys. B223 (1983), 218. 

[80] F. Bruckmann, T. Heinzl . T. Tok and A . Wipf, Instantons and Gribov copies in the 
maximally Abelian gauge, hep-th/0001175 . 



76 



(a) 













Fig. 9. Examples of Feynmann diagrams and the corresponding double line notations that appear 
in calculating the expectation value of the Wilson loop operator. The order of each diagram 
is estimated using the rule given in Fig. 8 as (a) Qym^^ = -^-^i (b) Sym^^ = A'^A'', (c) 
q^mN^ = A'A^, (d) q^mN = X'/N, (e) gf-M^V^ = A^iV, (f) gyM^V^ = A^iV, and (g) g^-uN'' = 
A^TV". [Here note that each contribution should be divided by the normalization factor A'^ for 
the fundamental quark, in agreement with the definition (3.6).] In the leading order of large 
A'' expansion, the leading contributions come from the planar diagrams, e.g., (a), (b), (c), (e) 
and (f), which are furthermore classified by the order of A. Note that the contribution from a 
nonplanar diagram (d) is suppressed for large A''. The diagram (g) is the vacuum polarization 
diagram due to quark-anti-quark pair creation and annihilation, which is neglected in the pure 
Yang-Mills theory without dynamical quarks. 
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(b) 



Fig. 10. A rectangular Wilson loop and the contribution to the Wilson integral in which x and y 
run over (a) opposite sides, and (b) same sides. 




Fig. 11. A circular Wilson loop. 





(c) 



A 



A 



(f ) 



Av 



A^ 




12. Graphical representation of the large expansion in the CP^ ^ model. The propagators 
are as follows: (a) (f> propagtor, Sabip^ + m^)~^; (b) a propagator, r{p)~^; (c) propagator, 

Sfiv — ^'^Ji" ^ [(p^ + 4:m^)f{p) — ^] ^ . The vertices are as follows: (d) a(j}a4>b vertex {-^Sab)', 

(e) Afj,(j)a(i>b vertex {^5ab{Pi^ +Pm)); (f) A^A^ffa'I'b vertex {^SabSf,,,)- 



79 



(a) (e) 




Fig. 13. Examples ol Feynmann diagrams, (a) and (b) are tadpole diagrams of order iV^/^ (c),(d) 
and (e) are vacuum polarization diagrams of order A'^". (f) and (g) are order N~^''^ diagrams. 



